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I.  THE  AUTONOMOUS  SUSTAINMENT  CONTAINER 

CONCEPT 


A.  BACKGROUND 

Throughout  the  history  of  sea  to  shore  engagements,  supplying  goods  from  units 
at  sea  to  forces  ashore  has  met  with  several  shortcomings.  Most  notably,  in  the  World 
Wars,  Korea,  and  the  recent  conflicts  in  the  Middle  East,  continuing  to  support  troops 
ashore  from  forces  afloat  carries  with  it  high  risk,  both  to  materiel  and  personnel.  The 
current  combat  environment  has  placed  further  strain  on  this  link  of  the  supply  chain  as 
non-national  forces  look  for  easier  and  faster  ways  of  harming  the  war  effort  than 
attacking  fortified  and/or  heavily  manned  emplacements  and  convoys.  A  new  method  for 
delivering  the  necessary  sustainment  goods  to  forces  ashore  without  endangering  lives  or 
providing  easy  targets  of  opportunity  for  these  opposing  forces  needs  be  found. 


B.  HISTORY  OF  SEA  TO  SHORE  SUSTAINMENT 

1.  Traditional  Methods 

To  provide  supplies  to  forces  on  shore,  the  method  used  throughout  history 
required  the  supply  ship  to  dock  with  a  physical  structure  designed  to  move  goods  from 
ship  to  shore.  As  ships,  especially  cargo  ships,  grew  in  size,  so  did  the  required  dock 
size.  The  dawn  of  the  20th  century  saw  cargo  ships  that  were  too  large  to  traditionally 
dock  that  moved  supplies  from  anchor  via  smaller  utility  boats.  This  method  still 
required  a  ship  to  berth  at  a  physical  dock  to  unload  supplies.  World  War  II  saw  the 
advent  of  supply  sustainment  without  the  physical  dock.  Army  ships  simply  pulled  up  to 
the  beach  and  offloaded  supplies.  This  method  of  beach  landing  to  resupply  troops  on  the 
beach  has  remained  the  preferred  mode  since  that  time.  As  seen  in  Figure  1,  U.S.  troops 
deploy  from  an  Army  LCV  on  D-Day  1944  compared  to  U.S.  and  Royal  troop  landing  in 
2007.  Although  the  commissioning  dates  of  these  two  vessels  are  significantly  different, 
the  methodology  remains  primarily  the  same.  The  same  is  true  about  the  delivery  of  parts 
and  supplies  from  sea  to  shore. 
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Figure  1.  U.S.  Marine  landing  at  Omaha  Beach,  1944  (left)  (From:  10)  Joint 

U.S.  and  Royal  Army  troop  landing  2007  (right)  (From:  11) 


2.  Non-traditional  Techniques 

There  are  points  in  history,  where  circumstances  prevented  the  use  of  tradition 
supply  delivery  methods.  Winds,  waves,  prevailing  currents  and  even  single  individuals 
were  used  to  circumvent  the  need  for  both  berth  and  beach,  or  even  for  the  ship  to  come 
in  contact  with  land  at  all.  “...the  Tokyo  Express  was  [having  difficulty  delivering 
supplies]... so  the  measures  now  being  taken  were  desperate  ones.  Destroyers  tried  to 
supply  these  Imperial  troops  by  making  high-speed  runs  to  Guadalcanal  and  dropping  off 
strings  of  lashed-together  drums  into  which  supplies  had  been  sealed.  Barges  from  the 
island  then  were  to  tow  these  drums  ashore.  The  troops  managed  to  retrieve  only  about  30 
per  cent  of  these  supplies  that  Tanaka's  destroyers  cast  upon  the  water.  Tanaka's 
destroyers  raced  down  The  Slot  on  3  December  and  dropped  strings  of  1,500  drums.  But 
the  island  troops  managed  to  haul  in  only  about  300  of  these  from  the  waters  off 
Tassafaronga  and  Segilau.  On  11  December,  Tanaka  himself  led  one  in  his  speedy 
Teruzuki...  and  dropped  some  1,200  drums  of  supplies  off  Cape  Esperance  and  then 
headed  north  again.  Only  some  of  the  drums  were  recovered  by  the  troops  ashore... ”[5] 
The  plight  of  the  Japanese  troops  continued  until  the  islands  were  finally  conquered  by 
U.S.  troops.  Further  stories  of  the  Tokyo  Express  include  using  divers  to  clandestinely 
recover  supplies  floated  by  Japanese  destroyers  and  midnight  barge  runs  through  U.S. 
patrolled  waters.  When  the  traditional  ways  begin  to  fail,  it  is  essential  to  make  every 
effort  to  solve  the  problem  with  something  new. 
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3.  Unmanned  Ships  in  History 

There  have  been  a  few  instances  of  unmanned  and  semi-autonomous  vessels 
taking  part  in  important  naval  operations  where  manned  vessels  simply  would  not  have 
been  as  capable.  One  of  the  most  famous  uses  of  unmanned  craft  occurred  during  the 
battle  between  the  Spanish  Armada  and  the  Queen’s  Fleet  off  of  Calais.  Sir  Francis 
Drake  and  John  Hawkins  hoped  to  pull  the  Armada  out  of  safe  harbor  using  the  then 
famous  hellbumers  of  Antwerp,  ships  that  had  been  made  into  tremendous  bombs,  as  the 
tool  in  which  to  bring  the  Spanish  fleet  into  combat.  “Drake  offered  a  ship  of  his  own, 
the  Thomas  of  Plymouth,  of  two  hundred  tons,  and  Hawkins  one  of  his;  and  with  the 
mounting  enthusiasm  six  more  were  recruited... the  captains  scattered  to  the  work  of 
getting  them  ready  and  stuffing  them  with  everything  available  to  make  a  blaze.”[3, 
p271]  “But  not  only  were  spars  and  sails  and  rigging  left,  since  the  ships  would  bear 
down  on  the  anchorage... under  full  sail,  but  all  the  ships’  guns  were  left  double-shotted 
to  go  off  when  the  fire  made  them  hot  enough,  to  add  to  the  terror  if  not  the  destruction  of 
the  enemy.”[3,  p271]  “...the  Spanish  anchorage  could  see  plainly  eight  tall  ships  with  all 
sails  set  and  lines  of  fire  beginning  to  run  up  their  rigging,  driving  straight  towards  them 
with  wind  and  tide.”[3,  p277]  “...the  line  was  driving  along  so  fast,  pushed  by  a  strong 
wind,  a  spring  tide  and  the  set  of  the  Channel  current  all  working  together... ”[3,  p277] 
“...the  double  shotted  guns  heated  now  white  hot,  began  to  go  off,  spraying  their  shot  at 
random  over  the  water  and  by  the  force  of  their  recoil  sending  up  a  fountain  of  sparks  to 
be  blown  among  the  boats... the  six  remaining  ships  swept  past  and  bore  down  on  the 
anchored  fleet,  the  sound  of  their  exploding  guns  heard  above  the  roar  of  the  flames  and 
the  fountains  of  sparks  shooting  skyward.”[3,  p277]  “  ...most  captains  simply  cut  their 
cables  and  ran  before  the  wind,  scattering,  some  here,  some  there,  as  if  they  were  as  much 
afraid  of  one  another  as  they  were  of  the  fire-ships.  The  strong  set  of  the  current  and  the 
rising  gale  swept  the  whole  disorderly  mob  out  through  the  straits  and  on  toward  the 
sands  of  the  Flemish  coast.”[3,  p278]  The  autonomous  fire-ships  did  little  physical 
damage  to  the  Spanish  fleet,  they  were  after  all  not  really  the  hellbumers  the  Spanish  had 
feared,  but  when  dawn  approached  the  Armada  was  in  complete  disarray.  Unable  to 
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reform,  they  were  harried  for  days  by  the  Queen’s  fleet  until  able  to  flee  northward.  This 
marked  the  end  of  the  Spanish  Order  and  set  about  the  beginning  of  the  English 
domination  of  the  sea. 

4.  The  Emphasis  on  Unmanned  Craft 

With  the  development  of  new  concepts  in  military  operations  and  reductions  in 
manpower  of  our  military  forces,  the  promotion  of  autonomous  systems  has  been  pushed 
to  the  forefront.  Recent  improvements  in  wireless  and  satellite  communications,  sensor 
and  data  relay  systems  and  remote  operability  has  made  unmanned  vehicles  an  available 
option  for  easing  the  manpower  burden  felt  by  much  of  the  military.  An  autonomous 
vehicle  is  a  self-piloted  vehicle  that  does  not  require  an  operator  to  navigate  and 
accomplish  its  tasks.  As  the  necessity  to  protect  human  life  grows,  the  need  for 
autonomous  vehicles  will  grow  as  well.  The  U.S.  Air  force  has  proven,  in  combat,  the 
range  of  uses  for  unmanned  vehicles.  The  Air  force  Predator  has  run  countless  sorties  in 
Afghanistan  and  Iraq  with  missions  ranging  from  simple  surveillance  to  offensive 
surgical  strikes.  The  pilot  does  all  this  from  a  remote  operating  station  well  away  from 
the  threat  region.  Building  off  of  the  advances  and  successes  of  the  Predator,  Global 
Hawk  will  fulfill  mission  roles  from  long  range  surveillance  and  tracking  to  bombing 
strikes  against  hardened,  buried  and  bunkered  targets.  Further  UAV  systems  are  in  the 
works  including  the  Boeing  stealth  UAV,  the  MQ-8B  Fire  Scout  helicopter  UAV,  and  the 
autonomous  programmable  Silver  Fox  and  Scan  Eagle  UAVs  developed  at  the  Naval 
Postgraduate  School.  These  advances  have  spread  to  surface  and  subsurface  craft  as 
well.  Current  systems  in  development  include  the  BAE  system  Protector  USV  for  Coast 
Guard  applications,  and  the  autonomous  Sea  Fox  USV  and  Remus  UUV. 

C.  THE  DIFFICULTIES  OF  THE  SEA-SHORE  INTERFACE 

The  sea  shore  interface  is  an  area  rife  with  hazards.  Shallow  and  shoal  waters  are 
home  to  shipping  dangers  such  as  reefs,  mine  hazard  areas,  capsizing  surf  conditions,  and 
underwater  rock  outcroppings.  The  shoals  are  further  threatened  during  combat  as  the 
threat  level  increases  with  enemy  activities.  For  the  Marines  at  Guadalcanal  the  problem 

was  service  support  facilities.  “Without  port  facilities,  supplies  reached  the  troops  only 
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after  a  series  of  time-consuming  and  labor-intensive  equipment  transfers.  Supplies  were 
first  unloaded  from  Navy  ships  offshore  into  lighters  for  the  trip  to  the  beach.  There 
American  service  support  personnel  transferred  the  tonnage  to  trucks  that  hauled  it  inland 
to  several  dumps  on  roads  under  construction.  From  the  dumps  supplies  had  to  be  hand 
carried,  by  both  Americans  and  native  laborers...  As  the  fighting  moved  further  inland 
the  distance  between  dumps  and  front  line  lengthened  and  road  building  could  not 
progress  as  fast  as  assault  units  advanced...  American  troops  temporarily  solved  the 
distribution  problem  by  using  the  many  streams  and  rivers  on  the  island.  Loading 
supplies  into  small  boats... Americans  pushed  the  craft  through  the  water  as  close  to  the 
tactical  units  as  possible.”[4,  p26-27]  The  1st  Marine  Division  faced  problems  from  both 
the  shoals  and  the  enemy  as  they  tried  to  get  supplies  onto  the  beach  at  Inchon.  “At 
Inchon,  the  lack  of  adequate  port  facilities  made  it  essential  that  the  naval  logistics  forces 
rapidly  move  supplies  ‘over  the  beach.”’[6,  p40]  While  the  UDT  cleared  obstacles  and 
the  Navy  Seabees  built  a  large  pontoon  dock  and  reconstructed  port  facilities,  eight  LSTs 
beached  themselves,  under  fire,  and  began  the  task  of  disgorging  tons  of  supplies  and 
munitions.  “After  the  evening  landing,  some  of  [Captain  Watson]  Singer’s  sailors  and 
[General  Henry]  Crowe’s  marines  had  to  unload,  organize  and  distribute  the  supplies 
from  the  eight  LSTs  temporarily  immobilized  on  Red  Beach... ”[6,  p40]  The  precarious 
condition  of  these  beached  vessels  can  be  seen  in  Figure  2,  but  the  surf  and  sand  wasn’t 
the  only  threat  to  these  vessels  or  their  crews.  “Enemy  snipers  in  Inchon  shot  at  Crowe’s 
leathernecks  and  Singer’s  sailors,  outlined  as  they  were  in  the  glare  of  the  floodlights 
installed  by  the  Seabees.  Nonetheless,  the  men  continued  to  toil  throughout  the  long 
night  to  accomplish  their  vital  mission.  All  of  the  LST  captains  withdrew  their  emptied 
ships  on  the  morning’s  high  tide.”[6,  p40] 
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Photo  U  8G-(j-42Q027  LSTs  unlcad  mg  at  Inchon.  15  September  1950 


Figure  2.  Four  “high-and  -dry”  LSTs  disgorge  supplies  onto  Red  Beach  on  or 

about  16  September  (1950)  (From:  6) 

D.  THE  ASCC  MISSION 
1.  Purpose 

Waterborne  logistics  delivery  systems  such  as  the  mechanized  Landing  Craft 

Utility  (LCU),  Landing  Craft  Air  Cushion  (LCAC),  and  High  Speed  Vessels  (HSV)  have 

been  used  for  logistical  delivery  for  many  decades.  All  three  of  those  types  of  delivery 

vehicles  are  very  capable  at  providing  logistic  support  with  speed  and  efficiency; 

however,  they  all  provide  the  opportunity  for  the  loss  of  human  life.  The  ASCC  is 

designed  to  complete  the  same  missions  with  equal  or  better  speed,  cost  and  efficiency,  to 

effectively  remove  the  need  to  task  those  particular  vessels  with  logistic  support 

missions.  The  overall  main  purpose  of  the  ASCC  is  to  replace  the  need  for  the  logistical 

interface  of  sea  basing.  With  the  reduction  of  sea  basing  comes  with  it  the  reduction  of 

bulk  cargo  constraints  and  cargo  breakouts,  manpower  and  the  addition  in  time  savings. 

The  mission  application  has  immense  appeal  to  Joint  Logistics  over  the  Shore  (JLOTS). 

The  ASCC  goes  further;  to  reduce  the  threat  to  crew  and  service  support  personnel 

aboard  both  the  sea  based  logistics  ship  and  the  delivery  vessel  itself.  The  reduction  in 

manpower  requirements,  enhancements  of  personnel  safety,  and  reduction  of  local 

infrastructure  requirements  are  its  greatest  attractors. 
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2.  Mission  Profile 

Figure  3  shows  a  hypothetical  coastline  where  troops  positioned  ashore  are 
expecting  logistic  support.  The  opposing  forces  have  mined  the  coastline  and  set  up 
surface  to  surface  missile  sights  overlooking  the  bay.  Under  conventional  sea-shore 
support  missions,  the  logistics  craft  would  have  to  close  to  within  a  few  miles  of  shore  to 
offload  the  supply  carriers.  These  supply  carriers,  LSVs,  LCACs,  etc.,  then  must 
navigate  through  the  mine  fields  as  seen  in  Figure  4.  This  scenario  would  put  a 
significant  number  of  lives  at  risk  and  could  easily  lead  to  the  loss  of  equipment,  vessels, 
and  personnel.  Using  the  ASCC,  the  logistics  vessel  remains  well  outside  of  weapons 
range  while  offloading  the  ASCCs.  The  container  would  then  travel  along  a  defined  route 
until  it  reaches  its  designated  offload  area.  Once  onshore,  the  ASCC  could  be  unloaded 
manually,  loaded  on  to  cargo  trucks  for  further  transport,  or  reloaded  with  retrograde  to 
return  to  the  logistic  ship  off  shore.  Figure  5  shows  the  ASCCs  route  bypassing  the  mine 
field  and  supplying  the  troops  on  the  beach  with  needed  supplies. 


Figure  3.  The  ASCC  Operating  Environment  (From:  14) 
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Figure  4.  The  Threat  to  Personnel  Using  Traditional  Methods  (From:  14) 


Figure  5.  The  ASCC  Delivery  System  (From:  14) 


3.  The  AEPCO  Design 

Advanced  Engineering  &  Planning  Corporation,  Inc.  (AEPCO)  has  proposed  a 
Autonomous  Sustainment  Cargo  Container  that  would  solve  the  aforementioned 
problems.  The  ASCC  aims  to  reduce  manpower  requirements  for  various  littoral  vessels 
as  well  as  enhance  personal  safety  by  transporting  cargo  containers  from  the  logistics 
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ship  to  shore  using  an  autonomous  system  attached  directly  to  the  container.  The  typical 
ISO  container  can  be  seen  in  Figure  6  and  comes  in  a  20  foot  and  40  foot  length.  Initial 
designs  have  used  to  20  foot  ISO  container  to  produce  the  vessel  in  Figure  7.  The 
concept  has  a  transit  speed  between  15  and  20  knots  and  enough  range  to  transport  goods 
from  the  logistics  ship  offshore,  follow  a  programmed  route  past  hazards,  and  arrive  for 
pick  up  by  troops  ashore.  AEPCO  is  currently  involved  in  concept  exploration  in  order  to 
answer  the  fuel/power  and  speed/range  questions. 


Figure  6.  Typical  ISO  Container  (From:  12) 


Figure  7.  AEPCO  Proposed  ASCC  system  (From:  13) 


E.  THE  PURPOSE  OF  THIS  ANALYSIS 

The  purpose  of  this  analysis  is  to  perform  a  proof  of  concept  on  the  AEPCO 
ASCC  design  with  an  emphasis  on  seakeeping.  The  ASCC  parameters  will  be  input  into 

a  computer  simulation/calculation  program  and  the  results  will  be  compared  with 
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indicated  ship  performance  parameters.  Based  on  the  results,  areas  of  concern  or  design 
shortfalls  will  be  further  scrutinized  with  an  eye  toward  making  design  improvement 
suggestions.  The  simulated  model  will  test  the  vessel  in  the  areas  of  pitch,  roll,  vertical 
acceleration,  lateral  acceleration,  mast  submergence,  propeller  racing  due  to  broaching, 
and  bow  slamming.  These  results  will  then  be  compared  to  probable  operational  areas  to 
determine  the  operational  limits  to  the  particular  operating  environment.  This  will  then 
allow  for  the  calculation  and  display  of  the  overall  operational  capability  of  this  vessel.  It 
is  not  the  intention  of  this  study  to  conduct  a  detailed  structural  analysis  of  the  ASCC  or 
create  a  finite  element  model  of  the  seas  and  their  interaction  with  the  vessel  structure. 
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II.  SEAKEEPING  ANALYSIS 


A.  INITIAL  SETUP 

1.  The  Proposed  Layout 

The  AEPCO  ASCC  is  designed  to  use  a  standard  20  foot  ISO  container.  It 
focuses  on  the  80%  loading  condition.  This  is  both  the  most  common  and  most 
economical  way  to  deliver  goods.  This  study  will  also  use  the  80%  as  the  basis  for 
analysis  but  will  also  add  the  50,  60,  and  70%  loading  conditions  for  comparison. 
Dimensions  for  the  20  foot  ISO  at  80%  load  can  be  seen  in  Table  1.  This  data  was  used 
to  build  the  spreadsheet  seen  in  Table  2.  The  ASCC  KG  calculation  table  is  designed  to 
solve  for  the  various  loading  conditions  based  on  a  high,  mid-box,  or  low  loaded  center 
of  gravity.  The  distance  from  the  keel  to  the  center  of  gravity  is  called  KG.  The  most 
common  loading  condition  is  the  mid-box  load  where  items  inside  the  container  are 
uniformly  distributed  and  the  center  of  gravity  remains  roughly  in  the  center  of  the  load. 
For  ease  of  calculations,  KG  values  were  rounded  to  the  nearest  10  for  input  into  the 
analysis  program.  For  the  engine  intake,  a  snorkel  system  has  been  designed  to  also  act 
as  a  mast  head  to  comply  with  Coast  Guard  Regulations.  The  mast  head  height  is 
determined  by  the  distance  from  which  the  mast  head  light  can  be  seen.  Coast  Guard 
regulations  require  that  partially  submerged,  unmanned  vessels  have  a  range  of  visibility 
of  three  miles. [8]  This  range  converts  to  a  required  minimum  mast  head  height  of  one 
meter.  The  one  meter  tall  snorkel  height  is  used  for  determining  the  mast  submergence 
event  occurrences. 


11 


Table  1.  Dimensions  of  a  20  ft  ISO  Container  at  80%  Load  (From:  13) 


1  20'  ISO  Under  Loading  Condition  1 

%  Loaded 

80 

A  [lbs] 

42316 

Length  (L)  [ft] 

19.875 

Beam  (B)  [ft] 

8 

Height  (H)  [ft] 

8 

A[LT] 

18.891 

▼  [ft3] 

661.188 

Draft  (T)  [ft] 

4.158 

Length-to- 

Beam  Ratio 

(L/B) 

2.48 

Beam-to-Draft 
Ratio  (B/T) 

1.92 

L'  [ft] 

35.875 

L'/B 

4.48 
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Table  2.  ASCC  KG  calculation 


Percent  Load 

High  Load 

0 

5 

10 

20 

30 

40 

50 

60 

70 

80 

Midbox  Load 

weight 

8002.78 

10147.36 

12291.93 

16581.09 

20870.24 

25159.39 

29448.55 

33737.70 

38026.85 

42316.00 

Low  Load 

weight  (LT) 

3.57 

4.53 

5.49 

7.40 

9.32 

11.23 

13.15 

15.06 

16.98 

18.89 

V 

125.04 

158.55 

192.06 

259.08 

326.10 

393.12 

460.13 

527.15 

594.17 

661.19 

T 

0.79 

1.00 

1.21 

1.63 

2.05 

2.47 

2.89 

3.32 

3.74 

4.16 

KB 

0.39 

0.50 

0.60 

0.81 

1.03 

1.24 

1.45 

1.66 

1.87 

2.08 

BM 

54.25 

42.79 

35.32 

26.19 

20.80 

17.26 

14.74 

12.87 

11.42 

10.26 

KG  (est) 

GM 

4.25 

50.40 

39.04 

31.68 

22.75 

17.58 

14.24 

11.94 

10.28 

9.04 

8.09 

4 

50.65 

39.29 

31.93 

23.00 

17.83 

14.49 

12.19 

10.53 

9.29 

8.34 

3.75 

50.90 

39.54 

32.18 

23.25 

18.08 

14.74 

12.44 

10.78 

9.54 

8.59 

3.5 

51.15 

39.79 

32.43 

23.50 

18.33 

14.99 

12.69 

11.03 

9.79 

8.84 

3.25 

51.40 

40.04 

32.68 

23.75 

18.58 

15.24 

12.94 

11.28 

10.04 

9.09 

3 

51.65 

40.29 

32.93 

24.00 

18.83 

15.49 

13.19 

11.53 

10.29 

9.34 

2.75 

51.90 

40.54 

33.18 

24.25 

19.08 

15.74 

13.44 

11.78 

10.54 

9.59 

2.5 

52.15 

40.79 

33.43 

24.50 

19.33 

15.99 

13.69 

12.03 

10.79 

9.84 

2.25 

52.40 

41.04 

33.68 

24.75 

19.58 

16.24 

13.94 

12.28 

11.04 

10.09 

2 

52.65 

41.29 

33.93 

25.00 

19.83 

16.49 

14.19 

12.53 

11.29 

10.34 

1.75 

52.90 

41.54 

34.18 

25.25 

20.08 

16.74 

14.44 

12.78 

11.54 

10.59 

1.5 

53.15 

41.79 

34.43 

25.50 

20.33 

16.99 

14.69 

13.03 

11.79 

10.84 

1.25 

53.40 

42.04 

34.68 

25.75 

20.58 

17.24 

14.94 

13.28 

12.04 

11.09 

1 

53.65 

42.29 

34.93 

26.00 

20.83 

17.49 

15.19 

13.53 

12.29 

11.34 

Metric 

Percent  Load 

0 

5 

10 

20 

30 

40 

50 

60 

70 

80 

weight 

8002.78 

10147.36 

12291.93 

16581.09 

20870.24 

25159.39 

29448.55 

33737.70 

38026.85 

42316.00 

weight  (LT) 

3.57 

4.53 

5.49 

7.40 

9.32 

11.23 

13.15 

15.06 

16.98 

18.89 

weight  (kg) 

3630.06 

4602.84 

5575.62 

7521.18 

9466.74 

11412.30 

13357.86 

15303.42 

17248.98 

19194.54 

V(nV*) 

3.54 

4.49 

5.44 

7.33 

9.23 

11.13 

13.02 

14.92 

16.81 

18.71 

T  (m) 

0.24 

0.30 

0.37 

0.50 

0.62 

0.75 

0.88 

1.01 

1.14 

1.27 

KB 

0.12 

0.15 

0.18 

0.25 

0.31 

0.38 

0.44 

0.50 

0.57 

0.63 

BM 

5.04 

3.98 

3.28 

2.43 

1.93 

1.60 

1.37 

1.20 

1.06 

0.95 

KG  (est) 

GM 

1.2954 

53.35 

41.99 

34.63 

25.70 

20.53 

17.20 

14.90 

13.23 

11.99 

11.04 

1.2192 

53.43 

42.07 

34.71 

25.78 

20.61 

17.27 

14.97 

13.31 

12.07 

11.12 

1.143 

53.50 

42.14 

34.78 

25.86 

20.69 

17.35 

15.05 

13.38 

12.14 

11.20 

1.0668 

53.58 

42.22 

34.86 

25.93 

20.76 

17.43 

15.12 

13.46 

12.22 

11.27 

0.9906 

53.66 

42.30 

34.94 

26.01 

20.84 

17.50 

15.20 

13.54 

12.30 

11.35 

0.9144 

53.73 

42.37 

35.01 

26.09 

20.91 

17.58 

15.28 

13.61 

12.37 

11.43 

0.8382 

53.81 

42.45 

35.09 

26.16 

20.99 

17.66 

15.35 

13.69 

12.45 

11.50 

0.762 

53.88 

42.52 

35.16 

26.24 

21.07 

17.73 

15.43 

13.76 

12.52 

11.58 

0.6858 

53.96 

42.60 

35.24 

26.31 

21.14 

17.81 

15.50 

13.84 

12.60 

11.65 

0.6096 

54.04 

42.68 

35.32 

26.39 

21.22 

17.88 

15.58 

13.92 

12.68 

11.73 

0.5334 

54.11 

42.75 

35.39 

26.47 

21.30 

17.96 

15.66 

13.99 

12.75 

11.81 

0.4572 

54.19 

42.83 

35.47 

26.54 

21.37 

18.04 

15.73 

14.07 

12.83 

11.88 

0.381 

54.27 

42.90 

35.54 

26.62 

21.45 

18.11 

15.81 

14.15 

12.91 

11.96 

0.3048 

54.34 

42.98 

35.62 

26.69 

21.52 

18.19 

15.89 

14.22 

12.98 

12.03 

2.  Modeling  Waves  and  Sea  States 

Vessel  motions  in  regular  waves  were  found  using  the  strip  theory  of  Salveson, 
Tuck,  and  Faltinsen  as  modified  by  Beck  to  include  the  surge  degree  of  freedom.  [3]  The 
linear  simultaneous  equation  that  must  be  solved  is 

7=1-6 

^  {co]{Mjk  +  Ajk)  +  io)eBjk  +  Cjk)qk  =  Fj  +F;D  where  £  1-6  are  surge,  sway,  heave,  roll, 

£=1,6 

pitch,  and  yaw  respectively;  M  is  the  mass  matrix,  A  is  the  added  mass  matrix,  B  is  the 
damping  matrix,  C  is  the  hydrostatic  restoring  force  matrix,  F1  is  the  Froude 
Krylovexciting  force  in  the  jth  mode  of  motion,  and  FD  is  the  diffraction  exciting  force  in 
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the  jth  mode  of  motion.  For  a  port  and  starboard  symmetric  ship  the  mass  matrix  is  given 
as 


M 


jk 


M 

0 

0 

0 

Mzc 

0 

0 

M 

0 

-Mzc 

0 

Mxc 

0 

0 

M 

0 

-Mxc 

0 

0 

-Mzc 

0 

^44 

0 

~he 

Mzc 

0 

-Mxc 

0 

hs 

0 

0 

Mxc 

0 

-he 

0 

he 

where  M  is  the  vessels  total  mass,  x,  y,  and  z  are  the  centers  of  mass  in  the  longitudinal, 
transverse  and  vertical  directions,  and  I  is  the  mass  moment  of  inertia  about  the 
subscripted  axis.  In  long  crested  random  seas  the  response  spectrum  is  determined  by  the 

linear  relationship  SQ  +(a>,/3)  =  \H(co,j3)f  Sj  +(co)  where  Si+(co)  is  the  one-sided 

incident  wave  spectrum,  H(co,(3)  is  the  response  amplitude  operator  for  waves  propagating 
in  direction  (3,  and  S0+(co,P)  is  the  one-sided  response  spectrum  due  to  waves  traveling  in 
direction  (i  A  “cosine- squared”  spreading  function  is  used  such  that  the  area  under  the 
spectrum  is  the  mean  square  value  of  the  process.  Several  different  wave  spectra  were 
used  by  the  modeling  program.  The  Pierson-Moskowitz  spectrum  as 

Sj  +(<»)  =  0.008  lg2  /  <»V”°'32(g///s®  } }  where  Hs  is  the  significant  wave  height,  ISSC 

spectrum  as  Sj  +  (&>)  =  O.WH2scoAm  /<yV”a44(<B,”/‘2') )  where  com  =2n  I  Tm  and 
Tm  =  mean  period ,  ITTC  spectrum  based  on  zero  crossings  as 

Sj  +  (<»)  =  (H2  /  47r)(a>ll  /  co5  )e(_(1/ ,a>) )  where  coz  =2 n!  Tz  and 

Tz  =  zero  crossing  period ,  ITTC  spectrum  based  on  modal  (or  peak)  period  as 

Sj  +(&>)=). 31 25H2co4p  /  co5)ei~l25(<apl<a)  ]  where  cop=2n!Tp  and 

Tp  =  period  of  spectral  peak ,  the  Ochi  6-parameter  spectrum  as 

Sj  +(co)  =  1/4 £  [{0.25(41,  +1)<}/1/ /T(Aj)]x[Hl  / atiUj  +i)]xe{“0-25(4^+x^/'i,)4} 

j~  1.2 

where  0Jmj  =  modal  frequency  of  jth  hump,  A.  =  shape  parameter  for  jth  hump  ,  and 
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Hsj  =  significant  wave  height  of  jth  hump ,  and 


JONSWAP  spectrum  as 


X  =  gx  I  Uf0 

Ul0  =  wind  speed  at  10  meters 
x  =  fetch  in  meters 


Sj  +  {co)  =  ag  co  e 


■5  f-1.25 (®„/ffl)4)  . 


'.ye 


,  a  =  0.076X"022 
where 

com=  2x3. 5X-03\g/Ul0) 


7  =  3.3 

<j  =  0.07  for  co  <  com 
<y  =  0.09  for  co>  com 


Further  refinement  of  the  wave  spectrum  is  done  using  the  defined  sea  state 
significant  wave  height  and  modal  wave  period  found  in  Table  3.  The  region  specific 
percent  probability  values  are  used  for  later  operability  calculations. 


Table  3.  Sea  State  Definition  and  Example  Region  Probabilities 


Sea 

Stale 

Significant  Wave 
Height  (m) 

Modal  Wave 
Period  (sec) 

Percent  Probability 
(Region  specific) 

1 

0.1 

6.3 

1 

2 

0.4 

6.5 

6 

3 

1.1 

7.5 

15 

4 

2.4 

8.8 

32 

5 

3.5 

9.7 

21 

6 

5.5 

12.4 

15 

7 

8.0 

15.0 

8 

8 

12.0 

16.5 

2 

3.  Standards  of  Measurement  for  Naval  Vessels 

The  ASCC  parametric  limits  were  based  on  a  combination  of  Naval  standards 
shown  in  Figure  8,  operational  requirements,  and  mechanical  limitations.  Due  to  the 
unique  nature  of  the  ASCC,  its  lack  of  crew,  radar,  and  sensitive  electronics,  motion 
extremes  such  as  pitch,  roll,  slam,  vertical  acceleration  and  acceleration,  and  lateral 
velocity  do  not  need  to  remain  as  fixed  as  the  values  listed  in  Figure  8.  Pitch,  roll, 
vertical  and  lateral  accelerations,  slam  are  generally  associated  with  crew  comfort.  These 
motions  contribute  to  seasickness.  For  the  ASCC,  seasickness  in  the  crew  is  not  an  issue 
and  the  parametrics  can  be  freely  modified  to  reflect  the  unique  nature  of  this  vessel. 
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These  malleable  limits  were  adjusted  as  follows:  pitch  limits  were  increased  to  10 
degrees,  roll  went  to  15  degrees,  slam  became  1.0  g’s,  vertical  acceleration  was  increased 
to  1.0  g’s,  vertical  velocity  went  to  3.0  m/s,  and  lateral  accelerations  became  0.5  g’s.  For 
propeller  racing  events  Figure  8  shows  a  limit  of  25/100  which  when  calculated  for  a 
design  trip  of  40  NM  leads  to  a  value  of  350  racing  events  per  hour.  A  new  category, 
mast  submergence  events,  was  added  due  to  the  damage  that  could  occur  should  the 
engine  intake  pull  a  slug  of  water  10  cm  or  greater  in  length.  For  this  study,  the  limit  was 
set  to  the  most  conservative  value,  one  event  per  trip. 
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Figure  8.  12  Prescribed  Values  for  Naval  Vessel  Performance  (From:  1) 

B.  DESIGN  ANALYSIS 
1.  SHIPMO 

The  program  used  to  model  the  waves  and  wave  responses  is  called 

SHIPMO. BM.  The  program  predicts  ship  motions  in  six  degrees  of  freedom  (6-DOF) 

and  five  components  of  the  shear  and  bending  moment  distributions.  It  is  based  on  the 

strip  theory  approach  of  Salvesen,  Tuck,  and  Faltinsen  which  was  extended  by  Beck  to 

include  the  surge  degree  of  freedom.  The  program  computes  both  regular  and  irregular 

waves  and  uses  a  cosine-squared  spreading  function  for  short  crested  seas.  The  program 

is  based  on  the  slender  body  theory,  which  assumes  a  long  slender  vessel  with  fine  ends. 
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The  version  used  for  this  analysis  has  been  updated  to  include  barge  type  forms  with 
blunt  ends  which  take  into  account  the  end  effects.  This  is  of  particular  note  for  the  box 
shaped  ASCC.  9,196  runs  were  made  using  loads  from  50  to  80  percent  in  10  percent 
intervals,  sea  states  from  1  to  8,  and  design  speeds  from  0  to  20  knots. 

a.  Data  Input 

Vessel  specific  data  and  run  dynamics  are  input  in  a  Fortran  77  format. 
The  ships  length,  breadth,  and  design  draft  in  meters  are  entered  in  at  line  3.  The  vessels 
hull  form  is  entered  as  a  series  of  discrete  points  as  seen  in  Figure  9.  KG  in  meters 
measured  from  the  waterline  is  entered  as  a  difference  with  draft  in  line  8.  The  range  of 
speeds  in  meters  per  second  is  entered  into  line  9.  The  range  of  headings  in  degrees  for 
seas  is  in  line  11.  A  sample  input  file  can  be  found  in  Appendix  B.  Every  change  in 
draft,  KG,  and  sea  state  becomes  a  new  input  file,  which  leads  to  a  large  number  of 
required  runs  to  cover  all  likely  combinations. 


6*  (2F 10 .4)  1VAXIS<I,IS>,  ZAX!S(|,IS),  1-1,  NT  I 

There  is  on«  dele  card  6*  for  each  input  point.  Start  at  keel  end  work  up. 


figure  A-T 
input  urrset  points 


The-  maximum  i?iDUt  ool  ms 

(NT)  is  IS  Per  stotlan 


Figure  9.  Input  Points  Format  for  SHIPMO  (From:  2) 

b.  Data  Extraction 

The  process  through  which  the  SHIPMO  program  calculates  the  modeling 
results  can  be  seen  in  the  Computer  Program  Flow  Chart  in  Figure  10  and  11.  Each  run 
in  SHIPMO  is  output  as  a  *.txt  table  and  *.emf  graphics  file.  These  outputs  are  imported 
and  collated  through  a  previously  created  Matlab  code  shown  in  Appendix  A.  This  code 

pulls  together  the  sea  state  and  load  results  and  builds  a  graphic  representation  as 
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compared  to  the  input  ship  limits  as  determined  by  ship  standards  in  part  A3  and  mission 
requirements.  The  results  come  out  as  a  polar  plot  with  ship  heading  toward  the  polar 
left,  so  waves  striking  the  ships  bow  will  be  seen  coming  from  the  270  degree  point  on 
the  polar  plot.  Speed  rings  are  measured  from  the  center  being  zero,  the  %  ring  at  5 
knots,  the  V2  ring  at  10  knots,  the  %  ring  at  15  knots,  and  the  outer  ring  at  the  max 
20knots.  Measured  results  that  exceed  the  design  parameters  are  noted  in  red,  the  various 
other  shades  note  degrees  to  which  the  vessel  is  approaching  those  limits. 


Computer  Program  Flow  Chart 


Figure  10.  The  SHIPMO  Computer  Program  Flow  Chart  (From:  2) 
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STOP 


Figure  11.  The  SHIPMO  Computer  Program  Flow  Chart  cont.  (From:  2) 

C.  RESULTS 

1.  Initial  Findings 

Preliminary  review  of  the  outputs  showed  a  startling  number  of  red  areas  in  the  1 5 
to  20  knot  speed  range  in  head  and  fore  quartering  seas.  Selected  results  can  be  found  in 
Appendix  B.  The  output  results  for  pitch,  roll,  slam,  mast  submergence,  propeller  racing, 
and  lateral  acceleration  were  review  for  potential  problems  in  the  design.  The  results 
were  also  put  together  to  help  determine  what  categories  were  of  more  particular 
importance  then  others.  Many  of  the  results  are  based  on  the  principles  of  ship  motion 
involving  crew  comfort  and  operational  windows  for  crew  laded  vessels.  Seasick  causing 
motions  like  lateral  and  vertical  accelerations  can  be  completely  ignored  due  to  the  lack 
of  crew  on  the  ASCC.  Motions  which  cause  seasickness  but  also  have  the  potential  to  do 
damage,  like  slamming,  can  also  be  shown  to  be  outside  the  scope  of  this  study  as  the 
structural  strength  or  placement  of  delicate  components  are  not  included.  Propeller 
racing  events  although  significant  are  less  of  a  concern  as  the  location  of  the  propeller 

was  not  known  and  the  data  was  estimated  for  a  propeller  20  and  30  cm  from  the  keel. 
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Additionally,  Caterpillar  diesels  are  extremely  hardened  against  this  type  of  marine 
environment  abuse.  After  careful  review,  pitch,  roll  and  mast  submergence  showed 
themselves  as  the  most  particularly  significant  results  obtained.  Mast  submergence  is 
especially  important  due  to  the  possible  damage  that  water  intrusion  can  do  to  the 
Caterpillar  diesel  [15].  Both  pitch  and  roll  motions  can  be  shown  to  be  the  root  cause  to 
both  propeller  racing  and  mast  submergence  events  and  are  the  driving  factors  to  the 
limits  to  the  ASCC’s  operational  capabilities. 

2.  Results  of  Particular  Importance 


a.  Mast  Submergence 

Based  on  documentation  provided  by  Caterpillar  Marine,  the  proposed 


supplier  of  the  ASCCs  prime  mover,  water  intrusion  of  any  amount  into  the  intake  would 
result  in  serious  damage  [15].  This  damage  is  not  relegated  to  simply  the  intake  but 
could  damage  the  manifold,  turbo,  exhaust  and  the  block  itself.  Additionally,  a  shortage 
of  combustion  air  due  to  a  water  slug  in  the  intake  tube  would  lead  to  loss  of  combustion 
air  and  a  stalling  of  the  engine.  Without  crew  onboard  the  ASCC  would  suffer  this 
damage,  stall,  and  be  unable  to  continue  its  mission.  This  study  took  this  particular  event 
into  account  for  analysis.  A  single  mast  submergence  event  on  a  50  NM  trip  was 
considered  unacceptable  and  can  be  seen  in  red  in  the  following  figures.  Bow  and  fore 
quarter  seas  produce  the  worst  results.  Sea  states  above  3  in  all  configurations  showed  at 
least  one  event  per  trip  at  transit  speeds  over  15  knots. 


Mast  Sub.  (/trip);  Loading  =  80  %;  KG  =  90  on;  SS  1  Range  =  1  NM;  Mast  =  100  cm 


0.8 


0.6 


0.4 


0.2 


180 


0 


Figure  12.  Mast  Submergence  at  80%  Load,  Sea  State  1,  Midbox  Load. 
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Figure  13. 


Figure  14. 


Figure  15. 


Mast  Submergence  at  80%  Load,  Sea  State  2,  Midbox  Load. 


Moot  Sub.  (/trip);  Loading  -  80  %;  KG  -  90  cn;  SS  3  Range  -  1  NM;  Moot  -  100  cm 


Mast  Submergence  at  80%  Load,  Sea  State  3,  Midbox  Load. 


Mast  Sub.  (/trip);  Loading  =  80  %;  KG  =  90  cm;  SS  4  Range  =  1  NM;  Mast  =  100  cm 


Mast  Submergence  at  80%  Load,  Sea  State  4,  Midbox  Load. 
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Figure  16. 


Figure  17. 


Figure  18. 


Mast  Submergence  at  80%  Load,  Sea  State  5,  Midbox  Load. 


uu  cm 

I 


Moot  Sub.  (/trip);  Loading  -  80  %;  KG  -  90  cm;  SS  G  Range  -  1  NM;  Mact  - 


071 

0.2 


Mast  Submergence  at  80%  Load,  Sea  State  6,  Midbox  Load. 


Mast  Sub.  (/trip);  Loading  =  80  %;  KG  =  90  cm;  SS  7  Range  =  1  NM;  Mast  =  100  cm 


Mast  Submergence  at  80%  Load,  Sea  State  7,  Midbox  Load. 
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Msst  Sub.  (/trip);  Loading  =  8D  %;  KG  =  90  cm;  SS  8  ;  Range  =  1  NM;  Mast  =  100  cm 


Figure  19.  Mast  Submergence  at  80%  Load,  Sea  State  8,  Midbox  Load. 

b.  Propeller  Racing 

As  previously  noted,  the  propeller  racing  events  are  noteworthy  simply  to 
show  the  coincidence  of  pitch  and  roll  to  overall  ships  operability.  As  with  mast 
submergence,  bow  and  fore  quartering  seas  produce  the  most  events  and  these  start  as 
low  as  sea  state  3. 


Figure  20.  Propeller  racing,  30  cm  above  keel  at  80%  Load,  Sea  State  1,  Midbox 

Load. 
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Figure  21. 


Figure  22. 


Figure  23. 


Propeller  Racing  (/hr);  Loading  =  8D  %;  KG  =  90  cm;  SS  2  ;  Prop.  =  30  cm 


Propeller  racing,  30  cm  above  keel  at  80%  Load,  Sea  State  2,  Midbox 

Load. 


Propeller  Racing  (/hr);  Loading  =  80  %;  KG  =  9C  cm;  SS  3  ;  Prep.  =  30  cm 


Propeller  racing,  30  cm  above  keel  at  80%  Load,  Sea  State  3,  Midbox 

Load. 


Propeller  Racing  (/hr);  Loading  =  80  %;  KG  =  9C  cm;  SS  4  ;  Prep.  =  30  cm 


Propeller  racing,  30  cm  above  keel  at  80%  Load,  Sea  State  4,  Midbox 

Load. 
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Figure  24. 


Figure  25. 


Figure  26. 


Propeller  Racing  (/hr);  Loading  =  80  %;  KG  =  90  cm;  SS  5  ;  Prop.  =  30  cm 


Propeller  racing,  30  cm  above  keel  at  80%  Load,  Sea  State  5,  Midbox 

Load. 


Propeller  racing,  30  cm  above  keel  at  80%  Load,  Sea  State  6,  Midbox 

Load. 


180 


210 


240 


300 


i 

I 


Prcpeller  Racing  (/hr);  Loading  =  80  %;  KG  =  90  cm;  SS  7  ;  Prop. 


150 


270 


250 

200 

150 

100 

50 


Propeller  racing,  30  cm  above  keel  at  80%  Load,  Sea  State  7,  Midbox 

Load. 
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Prcpeller  Racing  (/hr);  Loading  =  80  %;  KG  =  90  cm;  SS  8  ;  Prop.  =  30  cm 


Figure  27.  Propeller  racing,  30  cm  above  keel  at  80%  Load,  Sea  State  8,  Midbox 

Load. 


By  lowering  the  estimated  propeller  depth  by  10  cm,  the  following  figures 
are  output.  Although  not  readily  obvious  in  this  format,  the  change  made  had  little  effect 
on  the  likelihood  of  racing  events.  This  will  be  addressed  more  thoroughly  in  a  follow-on 
operability  study. 


Figure  28.  Propeller  racing,  20  cm  above  keel  at  80%  Load,  Sea  State  1,  Midbox 

Load. 
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Figure  29. 


Figure  30. 


Figure  3 1 . 


Propeller  racing,  20  cm  above  keel  at  80%  Load,  Sea  State  2,  Midbox 

Load. 


Propeller  Racing  (/hr);  Loading  =  80  %;  KG  =  9C  cm;  SS  3  ;  Prep.  =  20  cm 


Propeller  racing,  20  cm  above  keel  at  80%  Load,  Sea  State  3,  Midbox 

Load. 


Propeller  Racing  (/hr);  Loading  =  8D  %;  KG  =  90  cm;  SS  4  ;  Prop.  =  20  cm 


L 

•  '  250 
■  -200 

■  -150 

La 100 
m  so 


Propeller  racing,  20  cm  above  keel  at  80%  Load,  Sea  State  4,  Midbox 

Load. 
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Propeller  Racing  (/hr);  Loading  =  8D  %;  KG  =  90  cm;  SS  5  ;  Prop.  =  20  cm 


Figure  32.  Propeller  racing,  20  cm  above  keel  at  80%  Load,  Sea  State  5,  Midbox 

Load. 
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Figure  33. 


Figure  34. 


Prcpcllcr  Racing  (/hr);  Loading  -  80  %;  KG  ~  90  cm;  SS  6  ;  Prop.  -  20  cm 


Propeller  racing,  20  cm  above  keel  at  80%  Load,  Sea  State  6,  Midbox 

Load. 

Propeller  Racing  (/hr);  Loading  =  80  %;  KG  =  90  :m;  SS  7  ;  Prop.  =  2C  cm 


Propeller  racing,  20  cm  above  keel  at  80%  Load,  Sea  State  7,  Midbox 

Load. 
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Propeller  Racing  (/hr);  Loading  =  80  %;  KG  =  9C  cm;  SS  8  ;  Prep.  =  20  cm 


Figure  35.  Propeller  racing,  20  cm  above  keel  at  80%  Load,  Sea  State  8,  Midbox 

Load. 

c.  Pitch  and  Roll 

The  results  show  a  distinct  interaction  between  pitch  and  roll  and  mast 
submergence  events  and  propeller  racing  events.  Note  in  the  following  figures  that  both 
pitch  and  roll  degrees  go  in  extremis  with  waves  on  the  bow  or  the  fore  quarters.  It  is 
also  noteworthy  that  where  the  red  areas  of  the  pitch  and  roll  motions  overlap,  both  mast 
submergence  and  propeller  racing  go  red. 


Pitch  (deg)  Loading  =  80  %;  KG  =  90  cm;  SS  1 


Figure  36.  Pitch  at  80%  Load,  Sea  State  1,  Midbox  Load. 
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Figure  37. 


Figure  38. 


Figure  39. 


Pitch  (deg)  Loading  =  80  %;  KG  =  90  cm;  SS  2 


Pitch  at  80%  Load,  Sea  State  2,  Midbox  Load. 


Pitch  (deg);  Loading  =  80  %;  KG  =  90  cm;  SS  3 


Pitch  at  80%  Load,  Sea  State  3,  Midbox  Load. 


Pitch  (deg)  Loading  =  80  %;  KG  =  90  cm;  SS  4 


Pitch  at  80%  Load,  Sea  State  4,  Midbox  Load. 


31 


Figure  40. 


Figure  41. 


Figure  42. 


Pitch  (deg);  Loading  =  80  %;  KG  =  90  cm;  SS  5 


Pitch  at  80%  Load,  Sea  State  5,  Midbox  Load. 


Pitch  (deg);  Loading  =  80  %;  KG  =90  cm;  SS  6 


Pitch  at  80%  Load,  Sea  State  6,  Midbox  Load. 


=  80  %;  KG 
90  m 


Pitch  (deg);  Loading 


=  93  cm;  SS  7 


Pitch  at  80%  Load,  Sea  State  7,  Midbox  Load. 
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Figure  43. 


Figure  44. 


Figure  45. 


Pitch  (deg);  Loading  =  80  %;  KG  =  90  cm;  SS  8 


Pitch  at  80%  Load,  Sea  State  8,  Midbox  Load. 


Roll  (deg);  Loading  =  80  %;  KG  =  90  cm;  SS  1 


Roll  at  80%  Load,  Sea  State  1,  Midbox  Load. 


Roll  (deg);  Loading  =  80  %;  KG  =  90  cm;  SS  2 


Roll  at  80%  Load,  Sea  State  2,  Midbox  Load. 
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Figure  46. 


Figure  47. 


Figure  48. 


Roll  (deg);  Loading  =  80  %;  KG  =  90  cm;  S3  3 


Roll  at  80%  Load,  Sea  State  3,  Midbox  Load. 


Roll  (deg);  Loading 
120 


=  80  %;  KG  =  90  cm;  SS  4 


Roll  at  80%  Load,  Sea  State  4,  Midbox  Load. 


Roll  (deg);  Loading  =  80  %;  KG  =  90  cm;  SS  5 
30  10 

60 


210 


Roll  at  80%  Load,  Sea  State  5,  Midbox  Load. 
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Figure  49. 


Figure  50. 


Figure  51. 


Roll  (deg);  Loading  =  80  %;  KG  =  90  cm;  S3  6 


Roll  at  80%  Load,  Sea  State  6,  Midbox  Load. 


Roll  (deg);  Loading  =  80  %;  KG  =  90  cm;  SS  7 


Roll  at  80%  Load,  Sea  State  7,  Midbox  Load. 


Roll  (deg);  Loading  =  80  %;  KG  =  90  cm;  SS  8 


Roll  at  80%  Load,  Sea  State  8,  Midbox  Load. 
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D.  CONCLUSSIONS 


1.  Design/Mission  Shortcomings 

While  the  goal  of  the  ASCC  is  to  operate  in  sea  states  above  the  current  systems, 
Figures  52  and  53  indicate  that  this  may  not  be  entirely  possible.  By  putting  in  a  regional 
probability  of  occurrence  for  the  different  sea  states  and  overlaying  them  with  the 
operability  of  the  vessel  with  all  of  the  previously  mentioned  parametrics,  it  can  be  seen 
that  only  in  sea  states  1  and  2  is  there  a  low  occurrence  of  operational  interruptions.  At 
sea  state  3  the  ASCC  would  only  be  able  to  operate  less  than  65%  of  the  time.  At  sea 
state  4,  the  set  goal,  operability  drops  below  45%.  It  must  be  noted  that  these  figures  are 
for  all  speeds  and  headings.  It  would  be  possible  to  increase  the  operability  of  the  ASCC 
by  mandating  wave  heading  and  slower  speeds,  but  under  the  current  design  sea  state  4  is 
simply  too  damaging  to  be  able  to  successfully  operate.  Also  note  that  very  little 
difference  is  made  between  the  results  of  the  two  assumed  propeller  height  values. 


Loading  =  EO  %;  KG=  90  cm  Range  =  1  NM;  Mast  =  100  cm;  Prop.  =  30  cm 


Figure  52.  Overall  Operability  Index  for  80%  Midbox  Load  and  30  cm  propeller 

height. 
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Loading  =  £0  %;  KG=  90  cm  Range  =  1  NM;  Mast  =  100  cm;  Prop.  =  20  cm 


Figure  53.  Overall  Operability  Index  for  80%  Midbox  Load  and  20  cm  propeller 

height. 

2.  Improvement  Possibilities 

Although  the  predicted  response  is  deemed  satisfactory  for  the  intended  purpose 
of  the  ASCC,  it  seakeeping  performance  could  be  improved  as  follows:  (1)  Increased 
mast  height  in  order  to  reduce  mast  submergence  events.  Shown  below  are  mast 
submergence  events  for  mast  height  configurations  from  1  to  2  meters,  50%  loaded 
(return  load),  a  KG  of  50  cm  (upper  midbox),  and  at  the  design  goal  of  sea  state  4.  As 
can  bee  seen,  raising  the  mast  head  height  and  engine  intake  significantly  reduces  the 
number  of  submergence  events  down  to  very  rare  for  the  2  m  case.  This,  however,  might 
increase  the  center  of  gravity  of  the  system  and  would  need  further  testing.  (2)  Increased 
roll  damping  through  rudder  control.  This,  however,  might  increase  the  added  resistance 
due  to  rudder  action  and,  therefore,  reduce  the  available  payload.  (3)  Ballasted  outriggers 
in  order  to  decrease  roll  period.  This,  however,  might  increase  resistance  and  complicate 
launch  and  recovery,  and  limit  the  number  of  containers  that  can  be  carried  by  the 
logistics  ship. 
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Figure  54. 


Figure  55. 


Figure  56. 


Mast  submergence  at  the  initial  mast  height  of  1  m. 


Mast  Sub.  (/trip);  Loading  =  5D  %;  KG  =  50  cm;  SS  4  ;  Range  =  4  NVl;  Masl  =  120  cm 


Mast  submergence  with  a  mast  height  of  1.2  m. 


40  cm 

I 


Mast  Sub.  (/trip);  Loading  =  5D  %;  KG  = 


50  cm;  SS  4  ; 
90  in 


Range  =  4  NM;  Masl  = 


Mast  submergence  with  a  mast  height  of  1.4  m. 
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Figure  57. 


Figure  58. 


Figure  59. 


Mast  Sub.  (/trip);  Loading  =  50  %;  KG  =  50  cm;  SS  4  ;  Range  =  4  NM;  Mast  =  160  cm 


Mast  submergence  with  a  mast  height  of  1.6  m. 


Mast  Sub.  i/trip);  Loading  =  50  %;  <G  =  50  cm;  SS  4  ;  Range  =  4  NM;  Mast  =  180  cm 


Mast  submergence  with  a  mast  height  of  1.8  m. 


Mast  Sub.  (/trip);  Loading  =  5D  %;  KG  =  50  cm;  SS  4  ;  Range  =  4  NM;  Masl  =  200  cm 


Mast  submergence  with  a  mast  height  of  2  m. 
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3.  The  Next  Step 

More  detailed  seakeeping  calculations  could  be  performed  by  taking  into 
consideration  the  exact  shape  of  the  bow  and  stem  units,  along  with  experimental 
verification.  Detailed  graphical  modeling  would  also  be  helpful  in  further  determining 
root  causes  to  the  mast  submergence  and  propeller  racing  events.  Agencies  such  as  the 
David  Taylor  Institute  and  Naval  Surface  Warfare  Center  Carderock  could  be  engaged  to 
conduct  both  digital  FEM  and  CDF  modeling  and  physical  tank  model  experiments  to 
further  refine  these  results. 
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APPENDIX  A.  MATLAB  COLLATION  FUNCTION 


function  f cn_polar_speed_loop (Loading_Cond, KG, z_prop_keel, z_mast_top, range_NM) 
%  ASCC  Container  Data  Set 

%  Horizontal/Vertical  motions  -  Polar  plots 
%  Speed  in  knots  -  Heading  (180  head,  0  stern) 

%  Contour  plots  (heading/speed) 

%  Two  parameter  Bretchneider  spectrum  -  Long  crested  seas 
%  Dimensional  version  (Metric  units) 


%  General  Loop  -  meant  for  overnight  runs 
warning  off 
%  GENERAL  DATA 


lambda_min  =  5.000; 
lambda_max  =  505.000; 
delta_lambda  =  10.060; 
rho  =1024.000; 
zeta  =  1.000; 
L  =  6.058; 
g  =  9.810; 

%  Sea  state  information 


%  Min  wave  length  (m) 

%  Max  wave  length  (m) 

%  Wave  length  increment  (m) 

%  Water  density 
%  Regular  wave  height 

%  Reference  length  for  nondimensionalization 
%  Gravitational  constant 


sig_wave_height  =  [0.1  0.4  1.1  2.4  3.5  5.5  8.0  12.0]; 
modal_period_sec  =  [6.3  6.5  7.5  8.8  9.7  12.4  15.0  16.5]; 
percent_probab  =  [1  6  15  32  21  15  8  2] ; 

%  Set  limits  for  operations  ( PNA  Vol .  Ill) 

%  All  values  in  RMS  single  amplitudes 
%  at  Point-of-Interest  (POI)  where  applicable 


c_p_roll_max  =  15; 
c_p_pitch_max  =  10; 
c_p_ver_acc_max  =  1.0; 
c_p_lat_acc_max  =  0.5; 
c_p_bow_acc_max  =  1.0; 
c_p_bow_vel_max  =  3.0; 
c_p_Prop_max  =  350; 
c_p_mast_trip_max  =  1; 


%  Roll  (deg) 

%  Pitch  (deg) 

%  Abs .  Vert.  Accel,  at  amidships  (g's) 
%  Abs.  Lat .  Accel,  at  amidships  (g's) 

%  Slam  Accel,  at  bow  (g's) 

%  Relative  vert,  velocity  at  bow  (m/s) 
%  Propeller  racing  events  (per  hour) 

%  Mast  submergence  events  (per  trip) 


c _ p _ roll  = 

c_p_pitch  = 
c_p_ver_acc  = 
c_p_lat_acc  = 
c_p_bow_acc  = 
c_p_bow_vel  = 
c_p_prop  = 
c_p_mast_trip 


[0:0.5: c_p_roll_max] ; 

[0:1: c_p_pitch_max] ; 

[0:0.05: c_p_ver_acc_max] ; 
[0:0.02: c_p_lat_acc_max] ; 
[0:0.05: c_p_bow_acc_max] ; 
[0:0.10: c_p_bow_vel_max] ; 
[0:10: c_p_prop_max] ; 

=  [0:0.1: c_p_mast_trip_max] ; 


c_p_roll_max_string  = 
c_p_pitch_max_string  = 
c_p_ver_acc_max_string  = 
c_p_lat_acc_max_string  = 
c_p_bow_acc_max_string  = 
c_p_bow_vel_max_string  = 
c_p_prop_max_string  = 
c_p_ma  s  t_t r ip_max_s  t  r i ng 


num2str (c_p_roll_max) ; 
num2str (c_p_pitch_max) ; 
num2str (c_p_ver_acc_max) ; 
num2str (c_p_lat_acc_max) ; 
num2str (c_p_bow_acc_max) ; 
num2str (c_p_bow_vel_max) ; 
num2str (c_p_prop_max) ; 

=  num2str (c_p_mast_trip_max) ; 


%  Speed  range  (m/sec) 


V_min  =  0; 

V_max  =  10; 
delta_V  =  1; 

Speed_Range_string= ' All  Speeds ' ; 
Loading_Cond_string  =  num2str (Loading_Cond) ; 
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if  Loading_Cond  ==  80 
draft  =  1.27; 

end 

if  Loading_Cond  ==  70 
draft  =  1.14; 

end 

if  Loading_Cond  ==  60 
draft  =  1.01; 

end 

if  Loading_Cond  ==  50 
draft  =  0.88; 

end 

KG_string  =  num2str(KG); 

%  Point  for  motions  calculation 

z_point  =  0.00; 
x_point  =  0.00; 
y_point  =  0.00; 

z_prop_keel_string  =  num2str ( z_prop_keel ) ; 
z_prop  =  z_prop_keel/100  -  draft; 
z_prop_string  =  num2str ( z_prop) ; 
x_prop  =  -L/2; 
y_prop  =  0.00; 

z_mast  =  z_mast_top/100  +2.5  -  draft; 
z_mast_top_string  =  num2str ( z_mast_top) ; 
z_mast_string  =  num2str (z_mast) ; 
x_mast  =  -L/2; 
y_mast  =  0.00; 

range_string  =  num2str (range_NM) ; 


%  Loop  over  sea  state 


i_Sea_State=0  ; 

for  Sea_State=l : 1 : 8 , 

i_Sea_State=i_Sea_State+l ; 

[Loading_Cond  KG  z_prop_keel  z_mast_top  range_NM  Sea_State] 
S_S_string  =  num2str (Sea_State) ; 

HS=sig_wave_height (Sea_State) ; 

HS_string=num2str (HS) ; 

T_m=modal_period_sec (Sea_State) ; 

T_m_string  =  num2str (T_m) ; 
ome  g  a_m=  2  *  p i / T_m ; 

%  Set  up  file  reading  format. 


trigg  =  30; 
flloc  =  25; 
f 21oc  =  26; 
f 31oc  =  27; 
f 41oc  =  28; 
f 51oc  =  29; 
f 61oc  =  30; 


lambda 

lambdaoverL 

wavenumber 

omega 

wavenumber 
omega 
period 
f ilesize 
lambda  size 


lambda_min : delta_lambda : lambda_max; 
lambda/L; 

2 . 0*pi . /lambda; 
sqrt (wavenumber*g) ; 
wavenumber ' ; 
omega ' ; 

2 . 0*pi . / omega; 
size (lambda) ; 
trigg*f ilesize (2) ; 


beta_incr  =  10; 


d  Vector  of  wavelengths 

%  Wave  number 
%  Wave  frequency 
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iV=0 ; 

for  V  =  V_min : delta_V : V_max, 

iV  m  iV+1 ; 

V_string  =  num2str(V); 

V_calc  =  V;  %  m/sec 

V_knots  m  V/0.5144; 

time_hrs  =  range_NM/ (V_knots+0 . 01 ) ; 

ibeta  =  0; 

for  beta  =  0 : beta_incr : 3 60 ,  %  Loop  on  sea  direction 

if  beta  >  180 

beta=  360-beta; 
end 

beta_string  =  num2str (beta) ; 
beta_calc  =  beta*pi/180; 
ibeta  «  ibeta+1; 


%  Load  data  file  V  beta.txt 


load_f ilename  = 

strcat ( ' results \ ' ,  Loading_Cond_string, ' , KG_string, , V_string, 
filename  =  load (load  f ilename); 


omegae  =  omega-wavenumber*V_calc*cos (beta_calc) ; 

periode  =  2 . 0*pi ./omegae; 


%  Set  mass  matrix  elements 


Mil  =  filename (1 : trigg : lambda_size, 
M12  =  filename (1 : trigg : lambda_size, 
M13  =  filename (1 : trigg : lambda_size, 
M14  =  filename (1 : trigg : lambda_size, 
M15  =  filename ( 1 : trigg : lambda_size , 
M16  =  filename (1 : trigg : lambda_size, 
M21  =  filename (2 : trigg : lambda_size, 
M22  =  filename (2 : trigg : lambda_size, 
M23  =  filename (2 : trigg : lambda_size, 
M24  =  filename (2 : trigg : lambda_size, 
M25  =  filename (2 : trigg : lambda_size, 
M26  =  filename (2 : trigg : lambda_size, 
M31  =  filename (3 : trigg : lambda_size, 
M32  =  filename ( 3 : trigg : lambda_size , 
M33  =  filename (3 : trigg : lambda_size, 
M34  =  filename (3 : trigg : lambda_size, 
M35  =  filename ( 3 : trigg : lambda_size , 
M36  =  filename ( 3 : trigg : lambda_size , 
M41  =  filename (4 : trigg : lambda_size, 
M42  =  filename (4 : trigg : lambda_size, 
M43  =  filename ( 4 : trigg : lambda_size , 
M44  =  filename ( 4 : trigg : lambda_size , 
M45  =  filename (4 : trigg : lambda_size, 
M46  =  filename (4 : trigg : lambda_size, 
M51  =  filename ( 5 : trigg : lambda_size , 
M52  =  filename ( 5 : trigg : lambda_size , 
M53  =  filename (5 : trigg : lambda_size, 
M54  =  filename (5 : trigg : lambda_size, 
M55  =  filename ( 5 : trigg : lambda_size , 
M56  =  filename ( 5 : trigg : lambda_size , 
M61  =  filename ( 6 : trigg : lambda_size, 
M62  =  filename ( 6 : trigg : lambda_size, 
M63  =  filename ( 6 : trigg : lambda_size , 
M64  =  filename ( 6 : trigg : lambda_size , 
M65  =  filename ( 6 : trigg : lambda_size , 
M66  =  filename ( 6 : trigg : lambda_size, 


M  = 


[Mil 

M12 

M13 

Ml  4 

M15 

Ml  6 ; 

M2 1 

M2  2 

M2  3 

M2  4 

M2  5 

M2  6 ; 

M31 

M32 

M3  3 

M3  4 

M3  5 

M3  6 ; 

M41 

M42 

M4  3 

M4  4 

M4  5 

M4  6; 

M51 

M52 

M53 

M54 

M55 

M5  6 ; 

M61 

M62 

M63 

M64 

M65 

M66] 

1)  ; 
2)  ; 

3)  ; 

4)  ; 

5)  ; 

6)  ; 
1)  ; 
2)  ; 

3)  ; 

4)  ; 

5)  ; 

6)  ; 
1)  ; 
2)  ; 

3)  ; 

4)  ; 

5)  ; 

6)  ; 
1)  ; 
2)  ; 

3)  ; 

4)  ; 

5)  ; 

6)  ; 
1)  ; 
2)  ; 

3)  ; 

4)  ; 

5)  ; 

6)  ; 
1)  ; 
2)  ; 

3)  ; 

4)  ; 

5)  ; 

6)  ; 


_ *  f beta_string,  '  . txt ' ) 

Frequency  of  encounter 
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%  Added  mass  terms 


All  =  filename ( 6+1 : trigg : lambda_size , 1 ) ; 
A12  =  filename ( 6+1 : trigg : lambda_size, 2 ) ; 
A13  =  filename ( 6+1 : trigg : lambda_size, 3 ) ; 
A14  =  filename ( 6+1 : trigg : lambda_size ,  4 ) ; 
A15  =  filename ( 6+1 : trigg : lambda_size ,  5 ) ; 
A16  =  filename (6+1 : trigg: lambda_size, 6) ; 
A21  =  filename ( 6+2 : trigg : lambda_size, 1 ) ; 
A22  =  filename ( 6+2 : trigg : lambda_size ,  2 ) ; 
A23  =  filename ( 6+2 : trigg : lambda_size ,  3 ) ; 
A24  =  filename ( 6+2 : trigg : lambda_size, 4 ) ; 
A25  =  filename ( 6+2 : trigg : lambda_size, 5 ) ; 
A26  =  filename (6+2 : trigg: lambda_size, 6) ; 
A31  =  filename ( 6+3 : trigg : lambda_size , 1 ) ; 
A32  =  filename ( 6+3 : trigg : lambda_size ,  2 ) ; 
A33  =  filename ( 6+3 : trigg : lambda_size, 3 ) ; 
A34  =  filename ( 6+3 : trigg : lambda_size, 4 ) ; 
A35  =  filename ( 6+3 : trigg : lambda_size ,  5 ) ; 
A36  =  filename ( 6+3 : trigg : lambda_size ,  6 ) ; 
A41  =  filename ( 6+4 : trigg : lambda_size, 1 ) ; 
A42  =  filename ( 6+4 : trigg : lambda_size, 2 ) ; 
A43  =  filename ( 6+4 : trigg : lambda_size ,  3 ) ; 
A44  =  filename ( 6+4 : trigg : lambda_size ,  4 ) ; 
A45  =  filename ( 6+4 : trigg : lambda_size, 5 ) ; 
A46  =  filename (6+4 : trigg: lambda_size, 6) ; 
A51  =  filename ( 6+5 : trigg : lambda_size , 1 ) ; 
A52  =  filename ( 6+5 : trigg : lambda_size ,  2 ) ; 
A53  =  filename ( 6+5 : trigg : lambda_size, 3 ) ; 
A54  =  filename ( 6+5 : trigg : lambda_size, 4 ) ; 
A55  =  filename ( 6+5 : trigg : lambda_size ,  5 ) ; 
A56  =  filename ( 6+5 : trigg : lambda_size ,  6 ) ; 
A61  =  filename ( 6+6 : trigg : lambda_size, 1 ) ; 
A62  =  filename ( 6+6 : trigg : lambda_size, 2 ) ; 
A63  =  filename ( 6+6 : trigg : lambda_size ,  3 ) ; 
A64  =  filename ( 6+6 : trigg : lambda_size ,  4 ) ; 
A65  =  filename ( 6+6 : trigg : lambda_size ,  5 ) ; 
A66  =  filename (6+6 : trigg: lambda_size, 6); 


A 


[All  A12  A13  A14 
A21  A22  A23  A24 
A31  A32  A3 3  A34 
A41  A42  A43  A44 
A51  A52  A53  A54 
A61  A62  A63  A64 


A15  A16; . 
A25  A26; . 
A3 5  A3 6; . 
A4  5  A4  6 ;  . 
A55  A56; . 
A65  A66] ; 


%  Damping  terms 

Bll 

m 

filename 

(12  +  1 

B12 

= 

filename 

(12  +  1 

B13 

= 

filename 

(12  +  1 

B14 

m 

filename 

(12  +  1 

B15 

= 

filename 

(12  +  1 

B16 

= 

filename 

(12  +  1 

B21 

= 

filename 

(12+2 

B22 

= 

filename 

(12  +  2 

B23 

= 

filename 

(12+2 

B24 

= 

filename 

(12+2 

B25 

= 

filename 

(12+2 

B26 

= 

filename 

(12+2 

B31 

= 

filename 

(12  +  3 

B32 

= 

filename 

(12  +  3 

B33 

= 

filename 

(12  +  3 

B34 

= 

filename 

(12  +  3 

B35 

= 

filename 

(12  +  3 

B36 

= 

filename 

(12  +  3 

B41 

= 

filename 

(12  +  4 

B42 

= 

filename 

(12  +  4 

B43 

= 

filename 

(12  +  4 

B44 

= 

filename 

(12  +  4 

B45 

= 

filename 

(12  +  4 

trigg : lambda_size , 1 ) ; 
trigg : lambda_size, 2 ) ; 
trigg : lambda_size, 3 ) ; 
trigg : lambda_size ,  4 ) ; 
trigg : lambda_size, 5 ) ; 
trigg: lambda_size, 6) ; 
trigg : lambda_size, 1 ) ; 
trigg : lambda_size ,  2 ) ; 
trigg: lambda_size, 3) ; 
trigg : lambda_size, 4 ) ; 
trigg : lambda_size, 5 ) ; 
trigg: lambda_size, 6) ; 
trigg : lambda_size , 1 ) ; 
trigg : lambda_size ,  2 ) ; 
trigg : lambda_size, 3 ) ; 
trigg : lambda_size, 4 ) ; 
trigg : lambda_size, 5 ) ; 
trigg: lambda_size, 6) ; 
trigg : lambda_size, 1 ) ; 
trigg : lambda_size, 2 ) ; 
trigg: lambda_size, 3) ; 
trigg : lambda_size ,  4 ) ; 
trigg : lambda_size, 5 ) ; 
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B46  =  filename (12+4 : trigg : lambda_size, 6) ; 
B51  =  filename ( 12+5 : trigg : lambda_size, 1 ) ; 
B52  =  filename ( 12+5 : trigg : lambda_size ,  2 ) ; 
B53  =  filename (12+5 : trigg: lambda_size, 3) ; 
B54  =  filename (12+5 : trigg : lambda_size, 4 ) ; 
B55  =  filename ( 12+5 : trigg : lambda_size ,  5 ) ; 
B56  =  filename ( 12+5 : trigg : lambda_size ,  6 ) ; 
B61  =  filename ( 12+6 : trigg : lambda_size, 1 ) ; 
B62  =  filename (12+6 : trigg : lambda_size, 2 ) ; 
B63  =  filename ( 12+6 : trigg : lambda_size ,  3 ) ; 
B64  =  filename ( 12+6 : trigg : lambda_size ,  4 ) ; 
B65  =  filename (12+6 : trigg: lambda_size, 5) ; 
B66  =  filename (12+6 : trigg : lambda_size, 6) ; 

B  =  [Bll  B12  B13  B14  B15  B16; . . . 

B21  B22  B23  B24  B25  B26; . . . 

B31  B32  B33  B34  B35  B36; . . . 

B41  B42  B43  B44  B45  B46;... 

B51  B52  B53  B54  B55  B56; . . . 

B61  B62  B63  B64  B65  B66]  ; 

%  Hydrostatic  terms 

Cll  =  filename ( 18+1 : trigg : lambda_size, 1 ) ; 
C12  =  filename (18+1 : trigg : lambda_size, 2 ) ; 
C13  =  filename (18+1 : trigg: lambda_size, 3) ; 
C14  =  filename (18+1 : trigg : lambda_size, 4 ) ; 
C15  =  filename (18+1 : trigg : lambda_size, 5) ; 
C16  =  filename (18+1 : trigg : lambda_size, 6) ; 
C21  =  filename ( 18+2 : trigg : lambda_size, 1 ) ; 
C22  =  filename (18+2 : trigg : lambda_size, 2 ) ; 
C23  =  filename (18+2 : trigg : lambda_size, 3) ; 
C24  =  filename (18+2 : trigg : lambda_size, 4 ) ; 
C25  =  filename (18+2 : trigg: lambda_size, 5) ; 
C26  =  filename (18+2 : trigg : lambda_size, 6) ; 
C31  =  filename ( 18+3 : trigg : lambda_size, 1 ) ; 
C32  =  filename (18+3 : trigg : lambda_size, 2 ) ; 
C33  =  filename (18+3 : trigg: lambda_size, 3) ; 
C34  =  filename (18+3 : trigg : lambda_size, 4 ) ; 
C35  =  filename (18+3 : trigg : lambda_size, 5) ; 
C36  =  filename (18+3 : trigg : lambda_size, 6) ; 
C41  =  filename ( 18+4 : trigg : lambda_size, 1 ) ; 
C42  =  filename (18+4 : trigg : lambda_size, 2 ) ; 
C43  =  filename (18+4 : trigg : lambda_size, 3) ; 
C44  =  filename (18+4 : trigg : lambda_size, 4 ) ; 
C45  =  filename (18+4 : trigg: lambda_size, 5) ; 
C46  =  filename (18+4 : trigg : lambda_size, 6) ; 
C51  =  filename ( 18+5 : trigg : lambda_size, 1 ) ; 
C52  =  filename (18+5 : trigg : lambda_size, 2 ) ; 
C53  =  filename (18+5 : trigg: lambda_size, 3) ; 
C54  =  filename (18+5 : trigg : lambda_size, 4 ) ; 
C55  =  filename (18+5 : trigg : lambda_size, 5) ; 
C56  =  filename (18+5 : trigg : lambda_size, 6) ; 
C61  =  filename ( 18+6 : trigg : lambda_size, 1 ) ; 
C62  =  filename (18+6 : trigg : lambda_size, 2 ) ; 
C63  =  filename (18+6 : trigg : lambda_size, 3) ; 
C64  =  filename (18+6 : trigg : lambda_size, 4 ) ; 
C65  =  filename (18+6 : trigg: lambda_size, 5) ; 
C66  =  filename (18+6 : trigg : lambda_size, 6) ; 

C  =  [Cll  C12  C13  C14  C15  C16; . . . 

C21  C22  C23  C24  C25  C26; . . . 

C31  C32  C33  C34  C35  C36; . . . 

C41  C42  C43  C44  C45  C46;... 

C51  C52  C53  C54  C55  C56; . . . 

C61  C62  C63  C64  C65  C66] ; 


%  Total  exciting  forces 


Fl_amp  =  filename ( flloc : trigg : lambda_size, 5 ) 
F2_amp  =  filename (f21oc : trigg: lambda_size, 5) 
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F3_amp 
F4_amp 
F5_amp 
F  6_amp 


filename (f31oc:trigg: lambda_size, 5 ) ; 
filename (f 41oc : trigg : lambda_size, 5) ; 
filename ( f 51oc : trigg : lambda_size ,  5 ) ; 
filename (f61oc:trigg: lambda_size, 5 ) ; 


Fl_pha 

F2_pha 

F3_pha 

F4_pha 

F5_pha 

F6_pha 


filename ( f Hoc : trigg : lambda_size ,  6 ) ; 
filename (f21oc : trigg : lambda_size, 6) ; 
filename (f31oc:trigg: lambda_size, 6) ; 
filename (f41oc:trigg: lambda_size, 6) ; 
filename ( f 51oc : trigg : lambda_size ,  6 ) ; 
filename ( f 61oc : trigg : lambda_size ,  6 ) ; 


FI  =  Fl_amp . *exp (i*Fl_pha . *pi/180 . 0 ) ; 
F2  =  F2_amp . *exp (i*F2_pha . *pi/180 . 0 ) ; 
F3  =  F3_amp . *exp (i*F3_pha . *pi/180 . 0 ) ; 
F4  =  F4_amp . *exp (i*F4_pha . *pi/180 . 0 ) ; 
F5  =  F5_amp . *exp (i*F5_pha . *pi/180 . 0 ) ; 
F6  =  F6_amp . *exp (i*F6_pha . *pi/180 . 0 ) ; 


F  =  [F1;F2;F3;F4;F5;F6] ; 


%  Calculate  motions 


Allbar 

= 

- (omegae . 

-2)  .  * 

(Mil 

+ 

All) 

+ 

i*omegae . *B11 

+ 

Cll 

A13bar 

= 

- (omegae . 

"2)  .* 

(M13 

+ 

A13 ) 

+ 

i*omegae . *B13 

+ 

C13 

A15bar 

= 

- (omegae . 

"2)  .* 

(M15 

+ 

A15 ) 

+ 

i*omegae . *B15 

+ 

C15 

A22bar 

= 

- (omegae . 

-2)  .* 

(M2  2 

+ 

A22) 

+ 

i*omegae . *B22 

+ 

C22 

A24bar 

= 

- (omegae . 

-2)  .* 

(M2  4 

+ 

A2  4) 

+ 

i*omegae . *B24 

+ 

C24 

A2  6bar 

= 

- (omegae . 

"2)  .* 

(M2  6 

+ 

A2  6) 

+ 

i*omegae . *B2  6 

+ 

C26 

A31bar 

= 

- (omegae . 

"2)  .* 

(M31 

+ 

A31) 

+ 

i*omegae . *B31 

+ 

C31 

A33bar 

= 

- (omegae . 

-2)  .  * 

(M33 

+ 

A33) 

+ 

i*omegae . *B33 

+ 

C33 

A35bar 

= 

- (omegae . 

A2 )  .  * 

(M35 

+ 

A35) 

+ 

i*omegae . *B35 

+ 

C35 

A42bar 

= 

- (omegae . 

A2 )  .* 

(M42 

+ 

A42 ) 

+ 

i*omegae . *B42 

+ 

C42 

A44bar 

= 

- (omegae . 

A2 )  .* 

(M44 

+ 

A44) 

+ 

i*omegae . *B44 

+ 

C44 

A  4  6bar 

= 

- (omegae . 

A2 )  .* 

(M4  6 

+ 

A4  6) 

+ 

i*omegae . *B4  6 

+ 

C46 

A51bar 

= 

- (omegae . 

A2 )  .* 

(M51 

+ 

A51 ) 

+ 

i*omegae . *B51 

+ 

C51 

A53bar 

= 

- (omegae . 

A2 )  .  * 

(M53 

+ 

A53 ) 

+ 

i*omegae . *B53 

+ 

C53 

A55bar 

= 

- (omegae . 

A2 )  .* 

(M55 

+ 

A55 ) 

+ 

i*omegae . *B55 

+ 

C55 

A62bar 

= 

- (omegae . 

A2 )  .* 

(M62 

+ 

A62) 

+ 

i*omegae . *B62 

+ 

C62 

A64bar 

= 

- (omegae . 

A2 )  .* 

(M64 

+ 

A64) 

+ 

i^omegae . *B64 

+ 

C64 

A66bar 

= 

- (omegae . 

A2 )  .* 

(M66 

+ 

A66) 

+ 

i^omegae . *B66 

+ 

C66 

etal  c 

= 

(A33bar . * 

A55bar  -  A35bar. 

*A53bar )  ./.  .  . 

(Allbar . *A33bar . *A55bar  -  Allbar . *A35bar . *A53bar  -  A31bar . *A13bar . *A55bar  + 
A31bar . *A15bar . *A53bar  +  A51bar . *A13bar . *A35bar  -  A51bar . *A15bar . *A33bar ) . *F1  -  ... 

(A13bar . *A55bar  -  A15bar . *A53bar) ./ (Allbar . *A33bar . *A55bar  - 
Allbar . *A35bar . *A53bar  -  A31bar . *A13bar . *A55bar  +  A31bar . *A15bar . *A53bar  + 

A51bar . *A13bar . *A35bar  -  A51bar . *A15bar . *A33bar ) . *F3  +  ... 

(A13bar . *A35bar  -  A15bar . *A33bar) ./ (Allbar . *A33bar . *A55bar  - 
Allbar . *A35bar . *A53bar  -  A31bar . *A13bar . *A55bar  +  A31bar . *A15bar . *A53bar  + 

A51bar . *A13bar . *A35bar  -  A51bar . *A15bar . *A33bar ) . *F5 ; 

eta2_c  =  (A4 4bar . *A6 6bar  -  A4 6bar . *A64bar ) . / . . . 

(A22bar . *A4 4bar . *A6 6bar  -  A22bar . *A4 6bar . *A64bar  -  A42bar . *A2 4bar . *A6 6bar  + 
A42bar . *A2 6bar . *A64bar  +  A62bar . *A2 4bar . *A4 6bar  -  A62bar . *A2 6bar . *A44bar ) . *F2  -  ... 

(A24bar . *A66bar  -  A2 6bar . *A64bar ) . / (A22bar . *A44bar . *A66bar  - 
A22bar . *A4 6bar . *A64bar  -  A42bar . *A2 4bar . *A6 6bar  +  A42bar . *A2 6bar . *A64bar  + 

A62bar . *A2 4bar . *A4 6bar  -  A62bar . *A2 6bar . *A44bar ) . *F4  +  ... 

(A24bar .  *A4 6bar  -  A2 6bar . *A44bar ) . / (A22bar . *A44bar . *A66bar  - 
A22bar . *A4 6bar . *A64bar  -  A42bar . *A2 4bar . *A6 6bar  +  A42bar . *A2 6bar . *A64bar  + 

A62bar . *A2 4bar . *A4 6bar  -  A62bar . *A26bar . *A44bar) . *F6; 

eta3_c  =  -  (A31bar . *A55bar  -  A35bar . *A51bar ) . / . . . 

(Allbar . *A33bar . *A55bar  -  Allbar . *A35bar . *A53bar  -  A31bar . *A13bar . *A55bar  + 
A31bar . *A15bar . *A53bar  +  A51bar . *A13bar . *A35bar  -  A51bar . *A15bar . *A33bar ) . *F1  +  ... 

(Allbar . *A55bar  -  A15bar . *A51bar) ./ (Allbar . *A33bar . *A55bar  - 
Allbar . *A35bar . *A53bar  -  A31bar . *A13bar . *A55bar  +  A31bar . *A15bar . *A53bar  + 

A51bar . *A13bar . *A35bar  -  A51bar . *A15bar . *A33bar ) . *F3  -  ... 

(Allbar . *A35bar  -  A15bar . *A31bar) ./ (Allbar . *A33bar . *A55bar  - 
Allbar . *A35bar . *A53bar  -  A31bar . *A13bar . *A55bar  +  A31bar . *A15bar . *A53bar  + 

A51bar . *A13bar . *A35bar  -  A51bar . *A15bar . *A33bar ) . *F5 ; 

eta4_c  =  - (A42bar . *A6 6bar  -  A4 6bar . *A62bar ) . / . . . 
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A42bar . 

A22bar . 
A62bar . 

A22bar . 
A62bar . 


A31bar . 

Allbar. 
A51bar . 

Allbar . 
A51bar . 


A42bar . 

A22bar . 
A62bar . 

A22bar . 
A62bar . 


(A22bar . *A44bar . *A66bar  -  A22bar . *A4 6bar . *A64bar  -  A42bar . *A2 4bar . *A6 6bar  + 
*A2  6bar . *A64bar  +  A62bar . *A2 4bar . *A4 6bar  -  A62bar . *A2 6bar . *A44bar ) . *F2  +  ... 

(A22bar . *A6 6bar  -  A2 6bar . *A62bar ) . / (A22bar . *A44bar . *A66bar  - 
*A4 6bar . *A64bar  -  A42bar . *A2 4bar . *A6 6bar  +  A42bar . *A2 6bar . *A64bar  + 

*A24bar . *A4 6bar  -  A62bar . *A2 6bar . *A44bar ) . *F4  -  ... 

(A22bar . *A4 6bar  -  A2 6bar . *A42bar ) . / (A22bar . *A44bar . *A66bar  - 
*A4 6bar . *A64bar  -  A42bar . *A2 4bar . *A6 6bar  +  A42bar . *A2 6bar . *A64bar  + 

*A24bar . *A4 6bar  -  A62bar . *A26bar . *A44bar) . *F6; 
eta5_c  =  (A31bar . *A53bar  -  A33bar . *A51bar ) . / . . . 

(Allbar . *A33bar . *A55bar  -  Allbar . *A35bar . *A53bar  -  A31bar . *A13bar . *A55bar  + 
*A15bar . *A53bar  +  A51bar . *A13bar . *A35bar  -  A51bar . *A15bar . *A33bar ) . *F1  -  ... 

(Allbar . *A53bar  -  A13bar . *A51bar) ./ (Allbar . *A33bar . *A55bar  - 
*A35bar . *A53bar  -  A31bar . *A13bar . *A55bar  +  A31bar . *A15bar . *A53bar  + 

*A13bar . *A35bar  -  A51bar . *A15bar . *A33bar ) . *F3  +  ... 

(Allbar . *A33bar  -  A13bar . *A31bar) ./ (Allbar . *A33bar . *A55bar  - 
*A35bar . *A53bar  -  A31bar . *A13bar . *A55bar  +  A31bar . *A15bar . *A53bar  + 

*A13bar . *A35bar  -  A51bar . *A15bar . *A33bar ) . *F5 ; 
eta6_c  =  (A42bar . *A64bar  -  A44bar . *A62bar ) . / . . . 

(A22bar . *A4 4bar . *A6 6bar  -  A22bar . *A4 6bar . *A64bar  -  A42bar . *A2 4bar . *A6 6bar  + 
*A2 6bar . *A64bar  +  A62bar . *A2 4bar . *A4 6bar  -  A62bar . *A2 6bar . *A44bar ) . *F2  -  ... 

(A22bar . *A64bar  -  A24bar . *A62bar ) . / (A22bar . *A44bar . *A66bar  - 
*A4 6bar . *A64bar  -  A42bar . *A2 4bar . *A6 6bar  +  A42bar . *A2 6bar . *A64bar  + 

*A24bar . *A46bar  -  A62bar . *A2 6bar . *A44bar ) . *F4  +  ... 

(A22bar . *A4 4bar  -  A24bar . *A42bar ) . / (A22bar . *A44bar . *A66bar  - 
*A4  6bar  .  lk'A64bar  -  A42bar  .  *A2 4bar  .  *A6 6bar  +  A42bar  .  *A2  6bar  .  lk'A64bar  + 

*A24bar . *A4 6bar  -  A62bar . *A26bar . *A44bar) . *F6; 


%  Absolute  motions 


xil_point 

xi2_point 

xi3_point 

xi3_prop 

xi3_mast 

xi3  bow 


etal_c 
eta2_c 
eta3_c 
eta3_c 
eta3_c 
eta3  c 


+  eta5_c*z_point 
+  eta6_c*x_point 

-  eta5_c*x_point 

-  eta5_c*x_prop 

-  eta5_c*x_mast 

-  eta5_c,k'  (L/2  )  ; 


-  eta6_c*y_point ; 

-  eta4_c*z_point ; 
+  eta4_c,k'y_point  ; 
+  eta4_c*y_prop; 

+  eta4_c*y_mast ; 


%  Relative  motions/velocities/accelerations 


xi3 

bow  rel 

= 

xi3  bow  -  zeta; 

% 

relative 

vi3~ 

bow  rel 

= 

i*omegae . *xi3  bow  rel; 

% 

relative 

vi3^ 

point 

= 

i*omegae . *xi3  point; 

% 

absolute 

xi3 

point  rel 

= 

xi3  point  -  zeta; 

% 

relative 

vi3 

point  rel 

= 

i*omegae . *xi3  point 

_rel ; 

% 

relative 

ai3 

point 

= 

-i* (omegae . A2 ) . *xi3 

point; 

% 

absolute 

ai2 

point 

= 

-i* (omegae . A2 ) . *xi2 

point; 

% 

absolute 

vi3 

bow 

= 

i*omegae . *xi3  bow; 

% 

absolute 

ai3 

bow 

= 

-i* (omegae . A2 ) . *xi3 

bow; 

% 

absolute 

xi3 

prop  rel 

= 

xi3  prop  -  zeta; 

% 

relative 

xi3 

mast  rel 

= 

xi3  mast  -  zeta; 

% 

relative 

%  Begin  random  wave  calculations 

%  Long  crested  seas;  Bretschneider  spectral  formulation 


%  forward  at  POI 
%  lateral  at  POI 
%  vertical  at  POI 
%  vertical  at  prop 
%  vertical  at  mast 
%  vertical  at  bow 


vertical  motion  at  bow 
vertical  velocity  at  bow 
vertical  velocity  at  POI 
vertical  motion  at  POI 
vertical  velocity  at  POI 
vertical  accel  at  POI 
lateral  accel  at  POI 
vertical  velocity  at  bow 
vertical  accel  at  bow 
vertical  motion  at  prop 
vertical  motion  at  mast 


%  Sea  spectrum 

S  =  (1 . 25/4  )  *  (  (omega_m/N4  .  /omega  .  A5)  *  (HSA2  )  )  .  *exp  (-1 . 25*  (omega_mA4  .  /omega  .  A4  )  )  ; 

Se  =  S . /abs ( (1- (2 . 0/g) *omega*V_calc*cos (beta_calc) ) ) ;  %  Convert  S (w)  to  S(we) 

%  Define  response  spectra 

Seta3_c 
Seta4_c 
Seta5_c 
Sxil_point 
Sxi2_point 
Sxi3_point 
Sxi3_bow 
Sxi3_bow_rel 
Svi3_point 
Sai3_point 
Sai2_point 
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=  ( (abs (eta3_c) ) . A2) .*Se; 

=  ( (abs (eta4_c) ) . A2) .*Se; 

=  (  (abs  (eta5_c) )  . A2 )  . *Se; 

=  (  (abs (xil_point ) )  . A2 )  . *Se; 

=  ( (abs (xi2_point ) ) . A2 ) . *Se; 

=  ( (abs (xi3_point ) ) . A2 ) . *Se; 

=  ( (abs (xi3_bow) ) . A2 ) . *Se; 

=  (  (abs (xi3_bow_rel) )  . A2 )  . *Se; 
=  ( (abs (vi3_point ) ) . A2 ) . *Se; 

=  ( (abs (ai3_point ) ) . A2 ) . *Se; 

=  ( (abs (ai2_point ) ) . A2 ) . *Se; 


Svi3_bow 

Svi3_bow_rel 

Sai3_bow 

Svi3_point_rel 

Sxi3_point_rel 

Sxi3_point_rel_2 

Sxi3_prop_rel 

Sxi3_prop_rel_2 

Sxi3_mast_rel 

Sxi3  mast  rel  2 


( (abs (vi3_bow) ) .  A2 )  . *Se; 

( (abs (vi3_bow_rel) ) . A2) . *Se; 

( (abs (ai3_bow) ) .  A2) .  *Se; 

( (abs (vi3_point_rel ) ) .  A2 ) . *Se; 

( (abs (xi3_point_rel ) ) .  A2 ) . *Se; 

(omegae .  A2 ) .  *  ( (abs (xi3_point_rel ) ) .  A2 ) . *Se; 
( (abs (xi3_prop_rel ) ) .  A2 ) . *Se; 

(omegae .  A2 ) .  *  ( (abs (xi3_prop_rel ) ) .  A2 ) . *Se; 

( (abs (xi3_mast_rel ) ) .  A2 ) . *Se; 

(omegae .  A2 ) . * ( (abs (xi3_mast_rel ) ) .  A2 ) . *Se; 


%  Initializations 


Seta3_c_i  =  0; 
Seta4_c_i  =  0; 
Seta5_c_i  =  0; 
Sxil_point_i  =  0; 
Sxi2_point_i  =  0; 
Sxi3_point_i  =  0; 
Sxi3_bow_i  =  0; 
Sxi3_bow_rel_i  =  0; 
Svi3_point_i  =  0; 
Sai3_point_i  =  0; 
Sai2_point_i  =  0; 
Svi3_bow_i  =  0; 
Svi3_bow_rel_i  =  0; 
Sai3_bow_i  =  0; 
Svi3_point_rel_i  =  0; 
Sxi3_point_rel_i  =  0; 
Sxi3_point_rel_2_i  =  0; 
Sxi3_prop_rel_i  =  0; 
Sxi3_prop_rel_2_i  =  0; 
Sxi3_mast_rel_i  =  0; 
Sxi3  mast_rel  2  i  =  0; 


%  Integral  S(w)*|RAO|A2 


for  1=2 : 1 ; £ ilesize (2 ) , 

Seta3_c_i  =  Seta3 

*  abs ( (omegae ( I — 1 ) -omegae (I) ) ) ; 

Seta4_c_i  =  Seta4 

*  abs ( (omegae (1-1) -omegae (I) ) ) ; 

Seta5_c_i  =  Seta5 

*  abs ( (omegae ( I — 1 ) -omegae (I) ) ) ; 

Sxil_point_i  =  Sxil_ 

1) )  *  abs ( (omegae ( 1-1 ) -omegae ( I ) 

Sxi2_point_i  =  Sxi2_ 

1) )  *  abs ( (omegae (1-1) -omegae (I ) 

Sxi3_point_i  =  Sxi3_ 

1) )  *  abs ( (omegae ( 1-1 ) -omegae ( I ) 

Sxi3_bow_i  =  Sxi3_ 

*  abs ( (omegae ( I — 1 ) -omegae (I) ) ) ; 

Sxi3_bow_rel_i  =  Sxi3_ 

Sxi3_bow_rel ( 1-1 ) )  *  abs ( (omegae ( I 

Svi3_point_i  =  Svi3_ 

1) )  *  abs ( (omegae (1-1) -omegae (I ) 

Sai3_point_i  =  Sai3_ 

1) )  *  abs ( (omegae ( 1-1 ) -omegae ( I ) 

Sai2_point_i  =  Sai2_ 

1) )  *  abs ( (omegae (1-1) -omegae (I ) 

Sai3_bow_i  =  Sai3_ 

*  abs ( (omegae ( I — 1 ) -omegae (I) ) ) ; 

Svi3_bow_i  =  Svi3 

*  abs ( (omegae (1-1) -omegae (I) ) ) ; 

Svi3_bow_rel_i  =  Svi3 

Svi3_bow_rel (1-1) )  *  abs ( (omegae ( I 

Svi3_point_rel_i  =  Svi3 
Svi3_point_rel ( 1-1 ) )  *  abs ( (omegae ( I 

Sxi3_point_rel_i  =  Sxi3 
Sxi3_point_rel ( 1-1 ) )  *  abs ( (omegae ( I 


point_i 
’) ) ; 

point_i 
') )  ; 

point_i 
) ) ; 
bow_i 

bow_rel_i 
1 ) -omegae 
point_i 
') )  ; 

point_i 
’) ) ; 

point_i 
’) )  ; 
bow_i 

bow_i 

bow_rel_i 
■1)  -omegae 
point_rel 
■1 )  -omegae 
point_rel_: 
1) -omegae 


+  0, 

+  0. 

+  0, 

+  0. 

+  0, 

+  0. 

+  0, 

+  0. 
(I) ) ) ; 
+  0, 

+  0. 

+  0, 

+  0. 

+  0, 

+  0. 
(i)  )  )  ; 
i  +  0, 
(I) ) ) ; 
i  +  0. 
(i)  )  )  ; 


5* (Seta3_c (I) 

5* (Seta4_c (I) 

5* (Seta5_c (I) 

5* (Sxil_point (I) 

5* (Sxi2_point (I) 

5*  (Sxi3_point (I) 

5* (Sxi3_bow  (I) 

5* (Sxi3_bow_rel  (I) 

5* (Svi3_point (I) 

5* (Sai3_point (I) 

5* (Sai2_point (I) 

5* (Sai3_bow (I) 

5* (Svi3_bow (I) 

5* (Svi3_bow_rel (I) 

5* (Svi3_point_rel (I) 
5* (Sxi3_point_rel (I) 


+  Seta3_c  ( I  —  1 ) ) 

+  Seta4_c (1-1) ) 

+  Seta5_c ( I — 1 ) ) 

+  Sxil_point ( I- 
+  Sxi2_point ( I- 
+  Sxi3_point ( I- 
+  Sxi3_bow (1-1) ) 
+ 

+  Svi3_point ( I- 
+  Sai3_point ( I- 
+  Sai2_point ( I- 
+  Sai3_bow (1-1) ) 
+  Svi3_bow (1-1) ) 
+ 

+ 
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Sxi3_prop_rel_i  =  Sxi3_prop_rel_i  +  0 . 5* (Sxi3_prop_rel (I)  + 

Sxi3_prop_rel (1-1) )  *  abs ( (omegae ( 1-1 ) -omegae ( I ) ) ) ; 

Sxi3_mast_rel_i  =  Sxi3_mast_rel_i  +  0 . 5* ( Sxi3_mast_rel ( I )  + 

Sxi3_mast_rel ( 1-1 ) )  *  abs ( (omegae ( 1-1 ) -omegae ( I ))) ; 

end 


%  Second  moment  m_2 
for  1=2 : 1 : f ilesize (2 ) , 

Sxi3_point_rel_2_i  =  Sxi3_point_rel_2_i  +  0 . 5* ( Sxi3_point_rel_2 ( I )  + 
Sxi3_point_rel_2 (1-1) )  *  abs ( (omegae ( 1-1 ) -omegae ( I ))) ; 

Sxi3_prop_rel_2_i  =  Sxi3_prop_rel_2_i  +  0.5* ( Sxi3_prop_rel_2 (I)  + 

Sxi3_prop_rel_2 (1-1) )  *  abs ( (omegae (1-1) -omegae (I) ) ) ; 

Sxi3_mast_rel_2_i  =  Sxi3_mast_rel_2_i  +  0.5* ( Sxi3_mast_rel_2 (I)  + 

Sxi3_mast_rel_2 (1-1) )  *  abs ( (omegae ( 1-1 ) -omegae ( I ))) ; 

end 


%  RMS  values 


RMS_eta3_c 

RMS_eta4_c 

RMS_eta5_c 

RMS_xil_point 

RMS_xi2_point 

RMS_xi3_point 

RMS_xi3_bow 

RMS_vi3_point 

RMS_ai3_point 

RMS_ai2_point 

RMS_vi3_bow 

RMS_ai3_bow 

RMS_vi3_bow_rel 

RMS_xi3_bow_rel 

RMS_vi3_point_rel 

RMS_xi3_point_rel 

RMS_xi3_point_rel_2 

RMS_xi3_prop_rel 

RMS_xi3_prop_rel_2 

RMS_xi3_mast_rel 

RMS  xi3  mast  rel  2 


sqrt (abs (Seta3_c_i) ) ; 

57 . 3* sqrt (abs (Seta4_c_i) ) ; 

57 . 3* sqrt (abs (Seta5_c_i) ) ; 
sqrt (abs ( Sxil_point_i ) ) ; 
sqrt (abs ( Sxi2_point_i ) ) ; 
sqrt (abs ( Sxi3_point_i ) ) ; 
sqrt (abs (Sxi3_bow_i) ) ; 
sqrt (abs ( Svi3_point_i ) ) ; 
sqrt (abs ( Sai3_point_i ) ) ; 
sqrt (abs ( Sai2_point_i ) ) ; 
sqrt (abs (Svi3_bow_i) ) ; 
sqrt (abs (Sai3_bow_i) ) ; 
sqrt (abs ( Svi3_bow_rel_i ) ) ; 
sqrt (abs ( Sxi3_bow_rel_i ) ) ; 
sqrt (abs ( Svi3_point_rel_i ) ) ; 
sqrt (abs ( Sxi3_point_rel_i ) ) ; 
sqrt (abs ( Sxi3_point_rel_2_i ) ) ; 
sqrt (abs ( Sxi3_prop_rel_i ) ) ; 
sqrt (abs ( Sxi3_prop_rel_2_i ) ) ; 
sqrt (abs ( Sxi3_mast_rel_i ) ) ; 
sqrt (abs ( Sxi3_mast_rel_2_i ) ) ; 


RMS_eta3_c_vector (ibeta, iV) 
RMS_eta4_c_vector (ibeta, iV) 
RMS_eta5_c_vector (ibeta, iV) 
RMS_xil_point_vector (ibeta, iV) 
RMS_xi2_point_vector (ibeta, iV) 
RMS_xi3_point_vector (ibeta, iV) 
RMS_xi3_bow_vector (ibeta, iV) 
RMS_vi3_point_vector (ibeta, iV) 
RMS_ai3_point_vector (ibeta, iV) 
RMS_ai2_point_vector (ibeta, iV) 
RMS_vi3_bow_vector (ibeta, iV) 
RMS_vi3_bow_rel_vector (ibeta, iV) 
RMS_xi3_bow_rel_vector (ibeta, iV) 
RMS_ai3_bow_vector (ibeta, iV) 
RMS_vi3_point_rel_vector (ibeta, iV) 
RMS_xi3_point_rel_vector (ibeta, iV) 
RMS_xi3_point_rel_2_vector (ibeta, iV) 
RMS_xi3_prop_rel_vector (ibeta, iV) 
RMS_xi3_prop_rel_2_vector (ibeta, iV) 
RMS_xi3_mast_rel_vector (ibeta, iV) 
RMS_xi3_mast_rel_2_vector (ibeta, iV) 


RMS_eta3_c; 
RMS_eta4_c; 
RMS_eta5_c; 
RMS_xil_point ; 
RMS_xi2_point ; 
RMS_xi3_point  ; 
RMS_xi3_bow; 
RMS_vi3_point ; 
RMS_ai3_point ; 
RMS_ai2_point ; 
RMS_vi3_bow; 
RMS_vi3_bow_rel ; 
RMS_xi3_bow_rel; 
RMS_ai3_bow; 
RMS_vi3_point_rel ; 
RMS_xi3_point_rel ; 
RMS_xi3_point_rel_2 ; 
RMS_xi3_prop_rel; 
RMS_xi3_prop_rel_2 ; 
RMS_xi3_mast_rel ; 
RMS_xi3_mast_rel_2 ; 


N_P_prop (ibeta, iV) =3600* (1/sqrt (2*pi) ) * (RMS_xi3_prop_rel_2 /RMS_xi3_prop_rel ) *exp ( - 
2*z_propA2/ (2*RMS_xi3_prop_rel) ) ; 

N_P_mast (ibeta, iV) =time_hrs*3 600* (1/sqrt (2*pi) ) * (RMS_xi3_mast_rel_2 /RMS_xi3_mast_rel ) *exp 
(-2*z_mast A2/ (2*RMS_xi3_mast_rel) ) ; 
end 

end 
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lL=lambdaoverL ; 

temp_str  =  num2str (abs (str2num (beta_string) -180 ) ) ; 

%  Roll 
figure ( 1 ) 

[ th, r ] =meshgrid ( ( 0 : beta_incr : 3  6  0 ) *pi/180 , V_min : delta_V: V_max) ; 

[X, Y] =pol2cart (th, r) ; 
h=polar (th, r ) ; delete (h) ; 
hold  on 

contourf (X ' , Y ' , RMS_eta4_c_vector , c_p_roll) , colorbar 

title (['Roll  (deg);  Loading  =  '  Loading_Cond_string  '  %;  KG  =  '  KG_string  '  cm;  SS  ' 
S_S_string] ) 

f igure_string= [ ' graphs\roll\ ' , Loading_Cond_string, , KG_string, , S_S_string, ' . emf ' ] ; 

saveas (gcf , f igure_string, ' emf ' ) ; 

close 

%  Pitch 

figure (2) 

[ th, r ] =meshgrid ( ( 0 : beta_incr : 3  6  0 ) *pi/180 , V_min : delta_V : V_max) ; 

[X, Y] =pol2cart (th, r) ; 
h=polar (th, r ) ; delete (h) ; 
hold  on 

contourf (X ' , Y ' , RMS_eta5_c_vector , c_p_pitch) , colorbar 

title ([' Pitch  (deg);  Loading  =  '  Loading_Cond_string  '  %;  KG  =  '  KG_string  '  cm;  SS  ' 
S_S_string] ) 

f igure_string= [ ' graphs \pitch\ ' , Loading_Cond_string, ' _ 1 , KG_string, , S_S_string, 1 . emf ' ] ; 

saveas (gcf, f igure_string, ' emf ' ) ; 

close 

%  Absolute  Vertical  Acceleration  at  POI 
figure (3) 

[th, r] =meshgrid ( (0 :beta_incr : 360) *pi/180 , V_min : delta_V : V_max) ; 

[X, Y] =pol2cart (th, r) ; 
h=polar ( th, r ) ; delete (h) ; 
hold  on 

contourf (X ' , Y ' , RMS_ai3_point_vector/ g, c_p_ver_acc) , colorbar 

title ([' Absolute  Vertical  Acceleration  (g) ;  Loading  =  '  Loading_Cond_string  '  %;  KG  =  ' 
KG_string  '  cm;  SS  '  S_S_string] ) 

f igure_string= [ ' graphs \ver_acc\ ' , Loading_Cond_string,  ' _ ' , KG_string,  ' _ ' , S_S_string,  '  . emf  * ] 

saveas (gcf, f igure_string, ' emf ' ) ; 
close 

%  Absolute  Lateral  Acceleration  at  POI 
figure ( 4 ) 

[ th, r ] =meshgrid ( ( 0 : beta_incr : 3  6  0 ) *pi/180 , V_min : delta_V : V_max) ; 

[X, Y] =pol2cart (th, r) ; 
h=polar (th, r ) ; delete (h) ; 
hold  on 

contourf (X ' , Y ' , RMS_ai2_point_vector/g, c_p_lat_acc) , colorbar 

title ([' Absolute  Lateral  Acceleration  (g) ;  Loading  =  '  Loading_Cond_string  '  %;  KG  =  ' 
KG_string  '  cm;  SS  '  S_S_string] ) 

f igure_string= [ ' graphs \lat_acc\ ' , Loading_Cond_string, ' _ ' , KG_string, ' _ ' , S_S_string, ' . emf ' ] 

saveas (gcf, f igure_string, ' emf ' ) ; 
close 

%  Slam  Velocity  at  bow 
figure (5) 

[th, r] =meshgrid ( (0 :beta_incr : 360) *pi/180 , V_min : delta_V : V_max) ; 

[X, Y] =pol2cart (th, r) ; 
h=polar ( th, r ) ; delete (h) ; 
hold  on 

contourf (X ' , Y ' , RMS_vi3_bow_rel_vector , c_p_bow_vel ) , colorbar 
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title (['Slam  Velocity  (m/sec);  Loading  =  '  Loading_Cond_string  '  %;  KG  =  '  KG_string  ' 
cm;  SS  '  S_S_string] ) 

f igure_string= [ ' graphs \slam_vel\ ' ,  Loading_Cond_string, ' _ ' , KG_string, ' _ ' ,  S_S_string, ' . emf ' 

]  ; 

saveas (gcf ,  f igure_string, ' emf ' ) ; 
close 

%  Slam  Acceleration  at  bow 
figure ( 6 ) 

[th, r] =meshgrid ( (0 :beta_incr : 360) *pi/180 , V_min : delta_V : V_max) ; 

[X, Y] =pol2cart (th, r) ; 
h=polar ( th, r ) ; delete (h) ; 
hold  on 

contourf (X ' , Y ' , RMS_ai3_bow_vector/ g, c_p_bow_acc) , colorbar 

title ([' Slam  Acceleration  (g) ;  Loading  =  '  Loading_Cond_string  '  %;  KG  =  '  KG_string  ' 
cm;  SS  '  S_S_string] ) 

f igure_string= [ ' graphs\slam_acc\ ' , Loading_Cond_string, ' _ ' , KG_string, ' _ ' , S_S_string, '.emf' 

]  ; 

saveas (gcf, f igure_string, ' emf ' ) ; 
close 

%  Propeller  Racing  Events 
figure (7 ) 

[ th, r ] =meshgrid ( ( 0 : beta_incr : 3  6  0 ) *pi/180 , V_min : delta_V : V_max) ; 

[X, Y] =pol2cart (th, r) ; 
h=polar (th, r ) ; delete (h) ; 
hold  on 

contourf (X ' , Y ' , N_P_prop, c_p_prop) , colorbar 

title  ([' Propeller  Racing  (/hr);  Loading  =  '  Loading_Cond_string  '  %;  KG  =  '  KG_string  ' 
cm;  SS  '  S_S_string  '  ;  Prop.  =  '  z_prop_keel_string  '  cm']) 

f igure_string= [ ' graphs\race\ ' , Loading_Cond_string, ' _ ' , KG_string, ' _ ' , S_S_string, '_',z_prop 

_keel_string, ' .emf' ] ; 

saveas (gcf, f igure_string, ' emf ' ) ; 

close 

%  Mast  Submergence  Events 
figure (8) 

[th, r] =meshgrid ( (0 :beta_incr : 360) *pi/180 , V_min : delta_V : V_max) ; 

[X, Y] =pol2cart (th, r) ; 
h=polar ( th, r ) ; delete (h) ; 
hold  on 

contourf (X ' , Y ' , N_P_mast , c_p_mast_trip) , colorbar 

title (['Mast  Sub.  (/trip);  Loading  =  '  Loading_Cond_string  '  %;  KG  =  '  KG_string  '  cm;  SS 
'  S_S_string  '  ;  Range  =  '  range_string  '  NM;  Mast  =  '  z_mast_top_string  '  cm'  ] ) 

f igure_string= [ ' graphs\mast\ ' , Loading_Cond_string, ' _ ' , KG_string, ' _ ' , S_S_string, , range_ 
string, ' _ ' , z_mast_top_string, ' . emf ' ] ; 
saveas (gcf, f igure_string, ' emf ' ) ; 
close 

%  Operating  Envelopes  (for  area  calcs  and  01  processing) 

figure ( 9 ) 
colormap ( [ 0  0  0 ] ) 

[th, r] =meshgrid ( (0 :beta_incr : 360) *pi/180 , V_min : delta_V : V_max) ; 

[X, Y] =pol2cart (th, r) ; 
h=polar ( th, r ) ; delete (h) ; 
hold  on 

contourf (X ', Y ', RMS_eta4_c_vector, [ c_p_roll_max  c_p_roll_max] ) 
contourf (X' , Y' , RMS_eta5_c_vector , [ c_p_pitch_max  c_p_pitch_max] ) 
contourf (X' , Y' , RMS_ai3_point_vector/g, [ c_p_ver_acc_max  c_p_ver_acc_max] ) 
contourf (X' , Y' , RMS_ai2_point_vector/g, [ c_p_lat_acc_max  c_p_lat_acc_max] ) 
contourf (X' , Y' , RMS_vi3_bow_rel_vector , [ c_p_bow_vel_max  c_p_bow_vel_max] ) 
contourf (X' , Y' , RMS_ai3_bow_vector/g, [ c_p_bow_acc_max  c_p_bow_acc_max] ) 
contourf (X' , Y' ,N_P_prop, [ c_p_prop_max  c_p_prop_max] ) 
contourf (X ', Y ', N_P_mast, [ c_p_mast_trip_max  c_p_mast_trip_max] ) 

title (['01;  Loading  =  '  Loading_Cond_string  '  %;  KG  =  '  KG_string  '  cm;  SS  '  S_S_string 
';  Range  =  '  range_string  '  NM;  Mast  =  '  z_mast_top_string  '  cm;  Prop.  =  ' 
z_prop_keel_string  '  cm']) 
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f igure_string= [ ' graphs\OI\ ' ,  Loading_Cond_string, , KG_string, 1 , S_S_string, , range_st 

ring, ' _ ' , z_mast_top_string, ' _ ' , z_prop_keel_string, '.bmp']; 

saveas (gcf , f igure_string, ' bmp ' ) ; 

f igl_data=imread ( f igure_string) ; 

suml_data=sum (sum (f igl_data) ) ; 

hold  off 

y2_data=0 ; 

y3_data=l ; 

x2_dat a=4 190 4922; 

x3_data=5  9138077; 

x_data=suml_data ; 

a_data= (y3_data-y2_data) / (x3_data-x2_data) ; 

b_data= (y2_data*x3_data-y3_data*x2_data) / (x3_data-x2_data) ; 
y_da  t  a=a_da  t  a  *  x_da  t  a +b_da  t  a ; 
y_data=y_data+0 .01; 
if  y_data  >  1 
y_data=l ; 

end 

if  y_data  <  0 
y_data=0 ; 

end 

01=100*y_data; 

0I_string  =  strcat ( 1 01= ' , num2str (01 , 3 ) , ' % ' ) ; 

%text (-17 , 10 , 0I_string) 

0I_vector ( i_Sea_State ) =100*y_data; 

clear  figl_data 

close 

%  Sea  state  loop  ends 
end 

%  Compute  overall  operability  index 

001  =  100*sum ( (01_vector/100 ) . * (percent_probab/100 ) ) ; 

00I_string  =  strcat (' 001= num2str (001 , 3 )) ; 
figure ( 10 ) 

colormap([0  0  1 ] ) , bar (0I_vector ) , xlabel ( ' Sea  State ') , ylabel (' Operability  Index  (%) ' ) 
axis ([090  100] ) 

title ([' Loading  =  1  Loading_Cond_string  '  %;  KG  =  '  KG_string  '  cm;  Range  =  ’ 
range_string  '  NM;  Mast  =  '  z_mast_top_string  '  cm;  Prop.  =  '  z_prop_keel_string  '  cm’]) 
text ( 7 , 80 , 00I_string) 

f igure_string= [ ' graphs\00I\ ' , Loading_Cond_string, ' _ ' , KG_string, ' _ ' , range_string, , z_mas 

t_top_string, ' _ ' , z_prop_keel_string, ' .bmp']; 

saveas (gcf , f igure_string, ' bmp ' ) ; 

close 

clear  all 
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ASCC 


APPENDIX  B.  SHIPMO  SAMPLE  INPUT  FILE 


o  11 
0.0000 


o 

I 

o 

-0.8800 

0.2000\ 

-0.8800 

0.4000 

Syx0 . 8800 

0.6000 

-0V8800 

0.8000 

i 

o 

00 

1.0000  ' 

"^.8800 

1.2190 

-bv8800 

1.2190 

i 

o 

'O 

1.2190 

-0.6000 

1.2190 

-0.4000 

1.2190 

-0.2000 

1.2190 

0.0000 

-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 

-0.38 

9999.0 

1.0000 

000.1000 


Measurements  fore  and 
aft  of  amidships 


Number  of 
sections 


Table  of  offsets 


Increment  of  speeds 


Range  of  speeds 


-(Draft -KG) 


0.0000 

0.0  0.0 
05.0000  0505.0000 
180.1000  10.0000 


10.0000  00.0000  '  10.0000 


1 .0000 


Increment  of  headings 


Range  of  headings 
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APPENDIX  C 


NONCRITICAL  OUTPUT  RESULTS 


Figure  60.  Slam  acceleration  at  80%  midbox  load. 

| 

i 


Figure  61.  Slam  acceleration  at  50%  midbox  load. 
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Figure  62.  Lateral  acceleration  at  80%  midbox  load. 


Figure  63.  Lateral  acceleration  at  50%  midbox  load. 
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Figure  64. 


Figure  65. 


Vertical  acceleration  at  80%  midbox  load 


Vertical  acceleration  at  50%  midbox  load 
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Figure  66. 


Slam  velocity  at  80%  midbox  load 
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if  LDKh>)  =  id  v  KJS  ■  30  tmt  4.  Ring*  =  1  I'M.  Mrel  ■  itt  cm;  hop  6l.  Lutdng  "liiHSi#  cm;  ££  4.  Ihngt  =  1  8*4.  Mill  ■  iflfi  cm;  Prop  =  M  cm 
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